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SECTION  I 


INTRODUCTION 

The  skin  friction  and  heat  transfer  can  be  significantly  larger  for 
turbulent  flow  over  a  rough  surface  as  compared  with  an  equivalent  tur- 
bul'^nt  flow  over  a  smooth  sui^ace.  Many  surfaces  of  engineering 
Interest  are  rough  in  the  aerodynamic  or  hydraulic  sense  (to  be  defined 
in  the  following  sections).  Examples  of  systems  for  which  surface 
roughness  is  an  important  concern  are  re-entry  vehicles,  missiles, 
stores  carried  externally  on  high  performance  aircraft,  ships'  hulls, 
turbines,  heat  exchangers,  piping  networks,  and  atmospheric  flows.  In 
light  of  this  broad  Interest  and  applicability,  there  is  a  significant 
Interest  in  accurate  predictive  models  for  turbulent  flow  over  a  rough 
surface. 

1.  PREDICTION  APPROACHES 

Some  possible  schemes  for  predicting  the  momentum  and  ent'alpy 
fields  (and  by  consequence  the  skin  friction  and  heat  transfer)  for 
turbulent  flow  over  a  rough  surface  are: 

1.  Solve  the  complete  unsteady  Navler-Stokes  equations  in  a 
micro-grid. 

2.  Solve  the  time-averaged  Navier-Stokes  equations  in  a 
micro-grid. 

3.  Solve  the  time-averaged  Navier-Stokes  equations  in  a 
macro-grid. 

4.  Solve  the  time-averaged  boundary  layer  equations  in  a 
macro-grid. 

The  first  scheme  is  a  one-hundred  percent  first  principles  approach. 
The  complete  set  of  descriptive  equations  is  solved  numerically  in  a 
grid  system  that  is  fine  enough  to  resolve  the  details  of  both  the  tur¬ 
bulence  and  the  geometry  of  the  rough  surface.  Such  an  effort  would  be 
well  beyond  the  state  of  the  art.  The  computer  storage  and  computation¬ 
al  time  required  for  such  an  effort  would  be  astronomical. 

The  second  scheme  does  not  attempt  to  temporally  resolve  the  turbu¬ 
lent  fluctuations  but  does  use  a  grid  network  that  is  fine  enough  to 
resolve  the  geometry  of  the  rough  surface.  The  grid  for  such  a  scheme 
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would  still  be  beyond  the  present  computer  storage  capabilities.  For 
this  scheme  a  turbulent  closure  model  would  be  required  since  .he  time 
averaging  process  introduces  new  variables  (the  so-called  Reyno.'ds 
stresses)  but,  unfortunately,  no  new  equaLrons.  Therefore,  empirical 
i.nputs  are  required  to  close  the  system  of  equations. 

The  third  scheme  does  not  attempt  to  resolve  either  the  temporal 
details  of  the  turbulent  fluctuations  or  the  minute  details  of  the  rough 
surface  geometry.  That  is,  the  numerical  finite  difference  grid  is 
coarse  compared  with  the  scales  of  the  turbulence  and  of  the  roughness. 
This  relatively  course  grid  allows  the  problem  to  be  placed  on  current 
large  computers.  Here,  as  In  the  second  scheme,  a  turbulent  closure 
model  Is  required.  Also,  since  the  details  of  the  roughness  geometry  are 
not  resolved,  some  form  of  roughness  model  must  be  introduced  to  account 
for  the  effect  of  the  roughness  on  the  flow  field. 

The  fourth  set  is  actually  a  subset  of  the  third  and  is  mentioned 
here  for  completeness.  The  well-known  boundary  layer  approximations  lead 
to  a  simplified  version  of  the  Navler-Stokes  equations,  the  so-called 
boundary  layer  equations.  These  equations  are  parabolic  as  opposed  to 
the  elliptic  Navler-Stokes  equations  and  are  readily  solved  by  many  com¬ 
puter  codes.  Here,  as  with  the  third  scheme,  both  a  turbulence  closure 
model  and  a  roughness  model  are  required.  Because  of  the  slmpllf lcatio.< 
mentioned  above,  the  boundary  layer  equations  offer  comparative  computa¬ 
tional  efficiency.  Also,  many  flows  of  engineering  Interest  meet  the 
boundary  layer  assumptions.  Furthermore,  it  should  be  pointed  out  that 
as  turbulence,  closure  models  that  were  developed  for  boundary  layers  are 
readily  extended  to  turbulent  Navler-Stokes  calculations,  so  too  would  be 
roughness  models  that  are  developed  for  boundary  layers. 

From  the  above  discussion  it  is  apparent  that,  if  the  flow  over  a 
rough  surface  is  to  be  computed,  some  efficient,  accurate  roughness  model 
must  be  supplied.  There  are  two  main  schools  of  thought  on  the  modeling 
of  rough  surfaces  that  are  discusF.ed  in  detail  in  Section  II. 

They  are:  (1)  a  wholly  empirical  method  known  as  the  equivalent  sand-grain 
approach,  and  (2)  a  semi- empirical  method  known  as  the  discrete  element 
approach.  While  both  methods  require  experimental  inputs,  the  equivalent 
sand-grain  approach  may  require  experimental  data  on  the  particular 
surface  under  consideration.  On  the  other  hand,  the  discrete  element 
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approach  Incorporates  more  basic  physics  of  the  process  and  uses  a  more 
generalized  empirical  input.  It  is  therefore  applicable  to  a  broader 
spectrum  of  rough  surfaces  without  requiring  specific  experimental  data. 

2.  OBJECTIVE 

The  objective  of  the  present  work  was  to  develop  and  verify  a  pre¬ 
dictive  model  for  turbulent  flow  over  rough  surfaces  consisting  of  dis¬ 
tributed  (or  three-dimensional)  roughness  elements.  The  desired  model 
should : 

1.  Accurately  predict  both  skin  friction  and  heat  transfer. 

2.  Not  require  a  data  base  for  each  particular  surface  under 
consideration. 

3.  CONTENTS 

The  general  organization  of  the  work  presented  in  this  report  is 
described  below.  Section  II  is  a  discussion  of  the  background  of  predic¬ 
tive  methods  for  turbulent  flow  over  rough  surfaces.  The  present 
predictive  model  is  derived  and  discussed  in  Section  III.  In  Section  IV 
the  method  of  solution  is  presented  and  discussed.  The  calibration  of 
the  present  model  is  given  in  Section  V.  In  Section  VI  the  results  of 
the  present  model  are  compared  with  experimental  data  for  fully  developed 
flowr .  In  Section  VII  the  results  of  the  model  are  compared  with 
experimental  data  for  developing  boundary  layers.  The  conclusions  of  the 
work  and  recommendations  are  contained  in  Section  VIII. 

Ancillary  Information  is  contained  in  the  appendices. 
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SECTION  II 


BACKGROUND 


The  rigorous  study  of  the  effects  of  surface  roughness  on  fluid  flow 

and  heat  transfer  had  its  origin  with  the  classic  work  of  NikuraJse  [1] 

in  1933.  He  took  pressure  drop  and  velocity  profile  measurements  in 

pipes  roughened  with  tightly  sized  grains  of  "Goettingen"  sand.  He  made 

an  extensive  number  of  experimental  runs  c  werlng  6  sand  grain  size:: 

(r/k  ~  15,  30.6,  60,  126,  252,  507;  where  r  is  the  pipe  radius  and  k 
s  s 

is  the  sand  grain  size)  and  pipe  Reynolds  numbers  ranging  from  500  to 
1,000,000.  Inspection  of  Flg<<re  1,  which  shows  a  plot  of  Nlkuradse's 
friction  factor  versus  Reynolds  number,  reveals  the  magnitude  of  this 
data  base.  In  this  paper  Nlkuradse  Introduced  the  roughness  Reynolds 
number  as  a  measure  of  the  state  of  the  flow  over  the  rough  surface 
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where  u*  is  the  friction  velocity  (/x  /p)  and  k  is  the  size  of  the  sand 

W  8 

grains.  He  also  documented  the  three  regimes  of  rough  surface  flow: 


hydraulically  smooch  Re^^  <  5 

s 

transitionally  rough  5  <  RCj^  <55-70 

fully  rough  Re^  >55-70 

Hydraulically  (or  aerodynaolcally)  smooth  flow  is  flow  over  a  rough 
surface  chat  has  the  same  resistance  as  flow  over  a  smooth  surface  at  the 
same  Reynolds  number.  That  is,  the  surface  appears  smooch  to  the  fluid. 

It  is  characterized  by  the  skin  friction  coefficient  depending  on  the 
Reynolds  number  of  the  gross  flow  but  not  depending  on  the  roughness  size, 
shape,  density,  etc.  Fully  rough  flow  is  flow  where  Che  skin  friction 
coefficient  does  not  depend  on  the  Reynolds  number  of  the  gross  flow  and 
does  depend  solely  on  Che  character  of  the  roughness.  Transitionally 
'ough  flow  (net  Co  be  confused  with  the  transition  between  laminar  and 
turbulent  flow)  is,  as  its  name  suggests,  between  hydraulically  smooth 
and  fully  rough.  It  is  characterized  by  dependence  of  the  skin  friction 
coefficient  on  both  flow  Reynolds  number  and  roughness  character. 
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Figure  1.  Plot  of  Nlkuradse’s  (IJ  Friction  Factor  Data  Versus 
Pipe  Reynolds  Number 


These  three  regimes  are  usually  explained  in  terms  of  the  relation¬ 
ship  between  the  roughness  height  and  the  so-called  viscous  sublayer.  It 
is  said  that  in  the  hydraulically  smooth  regime  the  roughness  elements 
are  all  within  the  sublayer;  therefore,  viscous  effects  dominate  and  the 
surface  appears  smooth  to  the  flow.  In  the  transitionally  rough  regime 
the  roughness  elements  begin  to  protrude  through  the  sublayer  and  both 
viscous  and  roughness  effects  are  significant.  Finally,  in  the  fully 
rough  regime  the  sublayer  is  said  to  be  fully  destroyed  and  viscous 
effects  become  insignificant;  hence,  the  Reynolds  number  is  no  longer  a 
factor. 

The  sublayer  explanation,  while  enlightening,  depends  almost  com¬ 
pletely  on  the  height  of  the  roughness  and  ignores  the  Important  factors 
of  shape  and  spacing  density  of  the  roughness  elements.  An  otherwise 
smooth  surface  with  very  sparsely  spaced  large  (relative  to  the  sublayer) 
elements  may  appear  smooth  to  the  flow.  Morris  [2]  took,  a  different  view 
that  was  based  more  on  the  spacing  of  the  elements  than  the  height.  He 
classed  flow  over  a  rough  uurtace  into  four  regimes:  smooth,  semismooth, 
hyperturbulent  (wake  interference),  and  quasismooth.  The  smootn  regime 
is  equivalent  to  the  hydraulically  smooth  regime  discussed  above.  The 
semismooth  regime  consists  of  a  mostly  smooth  surface  with  very  sparsely" 
spaced  roughness  elements.  The  flow  resistance  is  viewed  as  a  sum  of 
the  viscous  resistance  on  the  smooth  parts  of  the  surface  and  the  pres¬ 
sure  drag  on  the  roughness  elements.  Hyperturbulent  flow  consists  of  a 
surface  where  the  elements  are  more  closely  spaced  and  the  wakes  behind 
the  elements  reach  the  adjacent  elements.  Quasismooth  flow  consists  of 
elements  so  closely  packed  that  the  crests  of  the  elements  and  the 
nearly  stagnant  fluid  in  the  intermediate  pockets  form  an  almost  smooth 
surface. 

Following  Nikuradse's  work,  the  next  major  step  in  the  development 
of  the  theory  of  flow  over  rough  surfaces  was  the  definitive  paper  by 
Schllchting  [3]  in  1936  Schlichtlng  conducted  experiments  in  a  rectan¬ 
gular  ciiannel  with  the  upper  surface  roughened  and  the  remaining  sides 
smooth.  His  purpose  was  to  study  the  effects  of  roughness  shape  and 
density  on  the  flow  resistance.  He  investigated  seven  different  shapes: 
large  spheres,  small  spheres,  spherical  segments,  cones,  short  angles, 
long  angles,  and  a  "Hamburg"  sand.  Each  shape  was  studied  over  a  range 
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of  spaclngs  and  a  limited  range  of  Reynolds  numbers.  Schllchtlng 's 
paper  is  definitive  for  two  reasons:  (1)  it  remains  the  most  comprehen¬ 
sive  set  of  controlled  experiments  for  flow  over  a  well-defined  rough 
surface,  and  (2)  two  methods  of  viewing  the  roughness  problem  were  put 
forward.  In  this  paper  Schlichting  developed  in  detail  the  well-known 
equivalent  sand-grain  roughness  approach  and  briefly  discussed  an  alter¬ 
native  approach  that  has  come  to  be  known  as  the  discrete  element 
approach. 

1.  EQUIVALENT  SAND-GRAIN  ROUGHNESS 

To  adequately  describe  a  rough  surface  at  three  measures  are 

required:  height,  spacing  density,  and  shape.  In  general,  a  multiplicity 

of  length  scales  are  required  (k, ,  k. ,  ...  k  ,  where  the  number  n  will 

1  /  n 

vary  from  surface  to  surface) .  Simplicity  has  been  a  powerful  incentive 

in  past  attempts  to  use  a  single  length  scale  to  describe  the  nature  of 

rough  surfaces.  The  parameter  that  is  most  commonly  used  is  the  well- 

known  equivalent  sand-grain  roughness,  k  . 

s 

In  his  paper  Schlichting  [3]  proposed  the  equivalent  sand-grain 

roughness  as  a  measure  of  the  flow  resistance  character  of  a  rough 

surface.  This  measure  was  a  length  scale,  k  ,  which  was  the  sand-grain 

s 

size  in  Nlkuradse's  experiment  that  gave  the  same  flow  resistance  at  the 
same  Reynolds  number  based  on  hydraulic  radius.  It  was  Schlichting 's 
stated  purpose  to  use  this  measure  only  as  a  means  of  extrapolating  a 
restricted  set  of  experimental  resistance  data  to  other  Reynolds  numbers 
based  on  Nlkuradse's  extensive  data  set.  To  quote  the  NACA  translation 
[3]: 

"It  is  the  object  of  our  test  to  be  able  to  specify  for 
each  type  of  roughness  Investigated  a  characteristic  number 
by  whose  aid  it  becomes  possible  to  predict  the  resistance 
for  the  same  type  of  roughness  at  other  Reynolds  numbers... 
than  those  directly  measured." 

"...It  appeared  advantageous  to  us  to  evaluate  our  results 
in  such  a  manner  that  they  could  easily  be  expressible  in  terms 
of  results  obtained  by  Nikuradse  for  sand  roughness  without, 
however,  necessarily  taking  the  latter  as  a  standard  type  of 
roughness.  One  objection  to  using  sand  roughness  as  a  standard 
is  that  this  type  is  not  satisfactorily  reproducible...." 

Schlichting 's  admonition  about  the  reproducibility  of  sand-type 

roughness  was  most  likely  prompted  by  his  experience  with  the  "Hamburg" 
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sand.  One  of  Schllchtlng's  roughness  media  was  a  sand  that  he  designated 
"Hamburg"  as  compared  to  the  "Goettingen"  sand  that  was  used  by 
Nikuradse.  The  "Hamburg"  sand  was  sized  in  the  same  manner  that 
Nikuradse  used  and  was  bonded  to  the  surface  in  the  same  manner  as 
Nikuradse.  For  the  "Hamburg"  sand  Schllchting  determined  an  equivalent 
sand  grain  roughness  k^/k  -  1.64  (where  here  k  is  the  size  of  the 
"Hamburg"  sand).  That  is,  the  "Hamburg"  sand  was  64%  rou'^ber  than  the 
same  size  "Goettingen"  sand.  Schllchting  also  tested  a  'Goettingen" 
sand  and  found  that  the  friction  factor  differed  (<nly  1  to  2%  from  those 
found  by  Nikuradse;  therefore,  the  difference  cannot  be  attributed  to 
differences  in  the  experimental  methods.  This  is  a  graphic  counterexample 
to  the  widely  held  belief  that  surfaces  chat  are  densely  roughened  with 
sand  or  that  appear  sand-graln-llke  (e.g.,  grit  blasted)  scale  on  a 
one-to-one  basis  with  Nikuradse 's  sand  roughness. 

The  equivalent  sand-grain  roughness,  k^,  is  usually  determined  by 
measuring  the  skin  friction  and  velocity  profiles  for  a  particular  sur¬ 
face  and  then  comparing  these  results  with  the  results  of  Nikuradse  [1]. 
For  example,  Schllchting  [3]  determined  his  values  for  k^  in  the  follow¬ 
ing  manner.  He  started  with  Che  logarithmic  velocity  profile  given  by 
Nikuradse  for  fully  rough  flow 

^  -  5.75  log  ^  +  8.48  (2) 

where  y  Is  the  distance  from  the  surface.  Then  he  assumed  chat  his 
velocity  profile  over  a  surface  of  elements  with  uniform  height,  k,  could 
be  fit  by  the  formula 

^  -  5.75  log  ^  +  A  (3) 

Schllchting  determined  the  friction  velocity  independently,  basically  by 
using  pressur .  drop  measurements.  The  intercept  A  was  determined  by 
graphically  fitting  equation  (3)  to  the  measured  velocity  profile. 

Then  by  setting  the  right-hand  sides  of  equations  (2)  and  (3)  equal 
to  each  other  he  obtained 

k 

5.75  log  -  8.48  -  A  (^a) 


or 
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k 

s 

k 


(8.48  -  A)/5.75 


(4b) 
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This  method  is  not  as  straightforward  as  it  appears.  The  reader  no 
doubt  recognizes  the  similarity  between  equations  (2)  and  (3)  and 
the  well-known  law  of  the  wall  for  turbulent  flow  over  smooth  surfaces 

5.75  legist  5.5  (5) 

This  similarity  is  not  coincidental.  In  their  formulations  of  equations 
(2)  and  (3)  Nikuradse  and  Schlichting  both  assumed  that  the  effects 
of  roughness  are  confined  to  a  region  very  near  the  wall;  therefore,  the 
slope  of  the  logarithmic  region  of  the  velocity  profile  is  the  same  for 
both  rough  and  smooth  surfaces  (the  only  difference  being  in  the  inter¬ 
cept).  This  argument  is  based  on  skin  friction  and  mean  velocity  mea¬ 
surements  along  with  a  .somewhat  nebulous  definition  of  the  origin  for 
the  y-axls.  Recent  extensive  measurements  (Pimenta  [4],  Coleman  [5]  and 
Ligranl  [6])  have  shown  that  roughness  affects  both  the  magnitude  and 
spectral  distribution  of  turbulent  kinetic  energy  across  the  entire 
boundary  layer,  suggesting  that  the  slope  of  the  logarithmic  region 
could  be  different  for  sioooth  and  rough  surfaces.  Furthermore,  Uram  [7] 
has  determined  logarithmic  region  slopes  ranging  from  4  to  9  for  differ¬ 
ent  rough  surfaces. 

The  appropriate  origin  for  the  distance  from  the  wall,  y,  must  be 
selected  if  the  velocity  profile  is  to  have  a  logarithmic  region  of  the 
form  of  equation  (3) .  Schlichting  used  the  melt  down  location  as 
the  origin  for  the  .  \xis  in  his  work.  That  is,  he  used  the  surface  that 
would  result  if  all  the  elements  were  melted  and  the  melt  was  allowed  to 
distribute  itself  evenly.  This  method  yields  a  hypothetical  surface  that 
maintains  the  same  fluid  volume.  This  location,  while  it  appears  logi¬ 
cal,  has  no  physical  basis  from  a  fluid  dynamics  viewpoint.  Figure  2 
shows  the  resulting  velocity  profile  plots  with  different  y-axls  origins 
for  a  surface  that  consists  of  densely  packed  spheres  [5].  Inspection 
of  the  figure  reveals  th«.t  different  (but  apparently  valid  from  visual 
inspection)  logarithmic  profiles  exist  for  each  of  the  y-axls  origins. 
Therefore,  the  mere  fact  that  a  logarithmic  region  seems  to  exist  does 
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Figure  2 . 


Typical  Velocity  Profiles  for  a  Fully  Rough  Turbulent 
Boundary  Layer  in  Zero  Pressure  Gradient  [5] 


not  always  suggest  that  the  proper  origin  has  been  selected.  A 
further  complication  is  the  fact  that  Nikuradse  did  not  state  the  origin 
for  y  in  his  report  [1].  Use  of  equation  (^)  requires  that  the  origins 
for  y  in  equations  (2)  and  (3)  be  the  same. 

Determination  of  the  correct  origin  for  y  is  one  of  the  more  diffi¬ 
cult  tasks  in  this  t3T>e  of  analysis.  Perry  and  Joubert  [8]  have  proposed 
a  graphical  method  for  determining  this  origin.  Also,  Monin  and  Yaglom 
[9]  have  proposed  a  method  based  on  the  mathematical  nature  of  the 
logarithm  function.  It  is  generally  believed  that  this  origin  should  lie 
somewhere  between  the  valleys  and  crests  of  the  elements.  However, 

Perry,  et  al  [10]  have  found  effective  wall  locations  1.5  to  2.0  element 
heights  below  the  crest  for  transverse  bar  roughness. 

If  the  equivalent  sand-grain  roughness  approach  is  to  be  a  useful 
predictive  trol,  some  method  must  exist  to  specify  k^  for  surfaces  for 
which  no  skin  friction  data  exist.  Dvorak  [11],  Simpson  [12],  and 

Dlrllng  [13]  have  presented  correlations  for  k  /k  (or  the  shift  in  the 

1  ® 

mean  velocity  profile,  Au  )  with  some  parameters  that  attempt  to  account 
for  the  dependence  on  roughness  height,  spacing,  and  shape.  The  correla¬ 
tion  of  Dlrllng  [13]  Is  representative  of  these  efforts  and  is  shown  In 

Figure  3.  Dlrllng  plotted  k  /k  versus  a  roughness  density  X  ■  (£  /k) 

1.25  ®  ^ 

(A^/A' )  ’  where  k  is  the  mean  roughness  height,  is  the  average 

center  to  center  spacing,  A  Is  the  windward  surface  area  of  the  element, 

8 

and  Ap  is  the  projected  windward  surface  area  of  the  element. 

Dlrllng 's  correlation  seems  to  collapse  the  available  data  better 

than  the  others.  However,  several  points  warrant  discussion.  Inspection 

of  Figure  3  reveals  that  the  data  for  k  /k  spread  significantly  for  a 

s 

given  value  of  For  X  •  25  it  is  seen  that  the  data  vary  from 

k  /k  5:  0.2  to  k  /k  =  0.6  and  for  X  ^  11  from  k  /k  s  0.8  to  k  /k  »  2.4, 
s  s  s  s 

while  the  correlation  gives,  respectively,  k  /k  »  0.3  and  k  /k  ••  1.2. 

s  s 

Therefore,  the  reliability  of  this  better  correlation  is  approximately 


^The  law  of  the  wall  for  rough  surfaces  is  often  expressed  alternatively 
as 
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Figure  3.  Effective  Sand-Grain  R'>*c‘'r.=-o  C..rrei 


±100%.  Berg  [14]  found  that  the  correlation  gave  a  value  of  three 

times  larger  than  his  experimentally  determined  values.  Secondly,  they 
do  not  distinguish  between  two-dimensional  and  three-dimensional  rough¬ 
ness.  The  nature  of  the  flow  around  these  two  different  types  of  rough¬ 
ness  is  very  different.  The  flow  over  two-dimensional  elements  is  more 
strongly  dependent  on  the  spacing  (distance  between  strips)  than  the 
height  while  for  three-dimensional  elements  both  height  and  spacing  seem 
to  be  about  equally  important.  For  example.  Perry,  et  al  [10]  found 
that  for  certain  ranges  of  k/2.^  the  skin  friction  coefficient,  C^,  does 
not  scale  with  the  roughness  height,  k,  for  transverse  bar  roughness. 
Also,  for  general  roughness  (i.e.,  randomly  sized  and  distributed 
elements)  there  is  no  set  correct  method  for  characterizing  the  roughness 
(i.e.,  specifying  k,  and  A^) .  King,  et  al  [15]  have  discussed 

this  problem  without  giving  a  resolution. 

All  of  the  above  correlations  rely  heavily  on  Schlichting's  [3] 

original  determinations  of  k  /k.  As  discussed  in  the  following  sections 

s 

and  in  Appendix  A,  Schlichting's  data  reduction  method  was  flawed 

(mainly  by  neglecting  the  effects  of  the  side  walls  of  his  4:1  aspect 

ratio  channel) .  As  shown  in  Appendix  A,  corrected  or  modified  values  of 

k  /k  differ  from  those  originally  repo*  ted  by  26%  to  555%. 
s 

The  equivalent  sand-grain  roughness,  k  ,  is  used  in  predictive 
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modeling  through  algebraic  correlations,  Integral  methods,  and  differen¬ 
tial  (finite  difference)  methods.  Many  correlatlona  for  skin  friction 
coefficient,  C^,  and  Stanton  number,  St,  have  been  presented.  Examples 
are  Schllchtlng  [3],  Dlpprey  and  Sabersky  [16],  Nestler  [17],  and  Seldman 
[18].  Dvorak  [11]  has  pr«>sented  an  integral  method  to  predict  skin 
friction  and  heat  transfer. 

In  the  present  state  of  boundary  layer  computation,  the  differential 
methods  are  by  far  the  most  important.  These  methods  solve  the  partial 
differential  equations  of  the  boundary  layer  with  finite  difference 
techniques.  Cebeci  and  Chang  [19]  present  a  method  that  relies  heavily 
on  the  equivalent  sand-grain  roughness.  They  solve  the  usual  incompress¬ 
ible  boundary  layer  continuity  and  momentum  equations  for  the  flow  field. 
They  account  for  the  roughness  via  a  virtual  origin  Ay  for  the  Prandtl 
mixing  length,  where  Ay  is  defined  as  a  function  of  k^.  One  of  the 
major  problems  with  differential  methods  that  use  the  equivalent 
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sand-grain  roughness  is  the  ill-defined  wall  boundary  condition.  Cebeci 
and  Chang  [19]  attempt  to  overcome  this  by  defining  wall  conditions  at 
some  distance  y^  above  the  crest  of  the  roughness.  They  determine  the 
velocity  boundary  conditions  this  point  from  empirical  velocity 

profile  correlations. 

Lin  and  Bywater  [20]  and  Adams  and  Hodge  [21]  present  differential 
methods  that  are  basically  discrete  element  approaches  (discussed  below) 
but  which,  however,  rely  on  equivalent  sand-grain  roughness  to  some 
degree  in  their  turbulence  models. 

The  reader  has  probably  noticed  that  the  above  discussion  of 
equivalent  sand-grain  roughness  is  almost  totally  concerned  with  momen¬ 
tum  transport.  In  fact, the  equivalent  sand-grain  roughness,  k^,  is 
defined  only  on  the  basis  of  skin  friction  and  velocity  profile  data. 

In  general,  attempts  to  use  equivalent  sand-grain  roughness  to  correlate 
heat  transfer  data  have  not  been  successful.  This  is  most  likely  because 
the  mechanisms  for  momentum  and  energy  transport  to  a  roughness  element 
are  different.  The  apparent  shear  stress  at  the  wall  (total  tangential 
force  on  the  wall  divided  by  wall  area)  is  composed  of  viscous  forces 
plus  form  drag  on  the  elements.  For  heat  transfer  there  is  no  transport 
mechanism  which  is  analogous  to  the  form  drag  on  the  element.  In  fact, 
there  is  no  physical  reason  for  two  surfaces  with  the  same  skin  friction 
coefficient  to  have  the  same  Stanton  number. 

In  suanary, the  outstanding  problems  with  the  equivalent  sand-grain 
approach  are  as  follows: 

1.  To  accurately  predict  skin  friction,  experimental  data  for 
the  particular  rough  surface  or  a  similar  surface  are 
usually  required. 

2.  In  general,  sand-like  surfaces  do  not  have  a  value  of 

k  /k  -  1. 

0 

3.  The  effective  wall  location  is  not  readily  apparent. 

4.  Many  roughness  geometries  do  not  correlate  with  a  single 
length  scale.  The  present  correlations  between  k^/k 
and  some  density  parameter  are,  in  general,  unreliable. 

5.  Attempts  to  Interface  equivalent  sand-grain  models  with 
finite  difference  boundary  layer  computations  results 
in  an  ill-defined  boundary  condition  at  the  wall. 
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6.  There  is  no  physical  reason  for  heat  transfer  data 

to  correlate  with  k  . 

s 

In  spite  of  appearances,  it  is  not  our  intention  here  to  say  that 
the  equivalent  sand-grain  approach  is  a  complete  failure.  Its  dominance 
of  the  subject  of  flow  over  rough  surfaces  for  the  past  half  century  is 
testimony  to  its  value  in  filling  some  of  the  gaps  in  the  science  of 
fluid  mechanics.  However,  many  gaps  remain.  In  recent  years  the 
alternate  approach  has  gained  momentum. 

2.  THE  DISCRETE  ELEMENT  APPROACH 

In  the  same  paper  [3]  in  which  Schlichtlng  proposed  the  equivalent 
sand-grain  roughness  approach,  he  proposed  (as  a  means  of  better  under¬ 
standing  the  dependence  of  flow  resistance  on  the  roughness  density) 
that  the  flow  resistance  of  the  rough  surface  be  divided  into  two  com¬ 
ponents:  (1)  that  due  to  the  form  drag  on  the  element  and  (2)  that  due  to 
the  viscous  shear  on  the  smooth  surface  area  between  the  roughness 
elements.  Schlichtlng  defined  W  as  the  total  resistance,  as  that 

part  of  W  due  to  the  drag  on  the  elements,  and  W  as  that  part  due  to  the 

8 

viscous  shear  on  the  smooth  parts  of  the  surface,  where  >  W  -  W  . 

^  2  ® 

Then  he  defined  an  element  resistance  coefficient  C.  •  2W  /(pu  ,  F  ), 

a  r  K  r 

where  u^^  is  the  velocity  at  the  crest  of  the  elements  and  is  the 
projected  frontal  area  of  the  roughness  element.  Schlichtlng  went  on  to 
calculate  C^  for  all  of  his  surfaces;  however,  this  Information  was  not 
put  to  any  practical  use  at  that  time.  Following  these  same  Ideas,  a 
body  of  literature  has.  arisen  that  here  is  referred  to  as  the  discrete 
element  approach. 

The  discrete  element  approach  considers  the  momentum  and  energy 
transport  processes  on  the  collection  of  individual  roughness  elements 
and  the  smooth  surface  between  the  elements.  The  basic  Idea  of  the 
discrete  element  approach  Is  to  formulate  a  system  of  partial  differen¬ 
tial  equations  that  describe  the  mass,  momentum,  and  energy  transport  for 
the  flow  over,  around,  and  between  the  roughness  "'r.  this  method 

the  roughness  effects  are  taken  as  an  integral  part  of  the  flow  problem 
and  not  (as  with  the  equivalent  sand  approach)  as  some  ill-defined 
boundary  condition. 
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■'igure  ^  shows  the  major  physical  phenomena  that  dlstingu  t.h  flow 
over  a  rough  surface.  They  are  blockage  of  the  flow,  form  drag  on  the 


roughness  elements,  and  local  heat  transfer  between  the  elements  and  the 
fluid.  The  roughness  elements  occupy  a  finite  fraction  of  the  tspace 
available  to  fluid  flow  and  thus  block  the  flow.  As  the  fluid  flows 
over  and  around  an  element,  zones  of  different  pressure  develop  resulting 
In  form  or  pressure  drag  on  the  elements.  Therefore,  because  of  these 
foices,  momentum  Is  extracted  volumetrlcally  from  the  flow.  If  the  fluid 
and  the  element  are  at  different  temperatures,  there  will  be  local  heat 
transfer  between  the  element  and  the  fluid.  Therefore,  the  elements  act 
also  as  volumetric  heat  sources  or  sinks  (depending,  of  course,  on  the 
relationship  between  the  fluid  and  element  temperatures). 

The  above  mentioned  partial  differential  equations  are  formulated 
in  such  a  way  that  they  account  for  the  roughness  geometry  and  the 
roughness  phenomena  (blockage,  distributed  momentum  sink,  and  distributed 
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heat  slnk/source) .  These  equations  are  then  solved  (usually  by  finite 
difference  methods)  vd.th  the  appropriate  boundary  conditions.  The  wall 
boundary  conditions  are  now  well  defined.  The  velocities  must  go  to 
zero  at  the  base  smooth  wall  and  the  temperature  of  the  fluid  at  the 
wall  must  be  equal  to  the  base  smooth  wall  temperature. 

The  other  common  wall  boundary  condition  for  the  energy  equation  is 
the  adiabatic  wall.  Care  must  be  taken  when  specifying  an  adiabatic 
rough  wall.  Merely  requiring  that  the  derivative  of  the  temperature  be 
zero  at  the  base  wall  is  not  adequate,  since  the  elements  would  remain 
as  a  major  route  for  energy  transport. 

The  discrete  element  approach  fills  many  of  the  gaps  left  by  the 
equivalent  sand-grain  approach.  A  brief  list  is  provided  below: 

1 .  The  roughness  geometry  is  accounted  for  in  the  formulation 
of  the  equations.  All  measures  are  included,  namely  height, 
shape,  and  density  or  distribution.  A  single  length  scale 

is  not  relied  upon.  Also,  no  nebulous  s..andard  roughness 
is  required. 

2 .  The  problem  of  the  effective  wall  location  and  the  ill- 
defined  boundary  conditions  are  not  encountered.  The 
discrete  element  partial  differential  equations  are 
Integrated  out  from  the  base  smooth  wall  while  the  rough¬ 
ness  effects  are  Included  in  the  equations. 

3 .  The  effect  of  roughness  on  energy  transport  is  included 
in  the  equations  in  a  physically  meaningful  manner. 

As  stated  abcve,  many  of  the  basic  ideas  of  the  discrete  element 
method  were  first  proposed  by  Schllchtlng  [3]  in  1936.  Later,  Llepmann 
and  Goddard  [22]  took  this  view  as  did  Lewis  [23].  In  recent  years 
several  serious  attempts  to  use  the  discrete  element  model  as  a  basis 
for  calculation  methods  have  been  presented.  Calculation  methods  for 
skin  friction  and/or  heat  transfer  on  rough  surfaces  using  the  discrete 
element  approach  have  been  reported  by  Adams  and  Hodge  [21],  Lin  and 
Bywater  [20],  Flnson  and  co-workers  [24,25,26],  Christoph  and  co-workers 
[27,28],  and  Roberson  and  co-workers  [29,30,31]. 

The  model  of  Adams  and  Hodge  [21]  includes  only  one  aspect  of  the 
discrete  element  approach,  namely  the  form  drag  of  the  roughness  elements. 
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They  do  not  address  the  problem  of  flow  blockage  caused  by  the  elements. 
Also,  they  do  not  model  the  heat  transfer  between  the  elements  and  the 
fluid.  However,  they  do  calculate  heat  fluxes  and  by  consequence 
Stanton  numbers  with  apparent  success  for  some  cases.  In  addition,  they 
only  use  the  discrete  element  approach  to  calculate  the  flow  field  over 
fully  rough  surfaces.  They  calculate  the  flow  field  over  transitionally 
rough  surfaces  by  using  an  equivalent  sand-grain  approach  from  Healzer 
[32]. 

Lin  and  Bywater  [20]  address  all  three  roughness  phenomena  (blockage, 
distributed  momentum  sink,  and  distributed  heat  sink/source).  However, 
they  do  not  distinguish  between  flow  blockage  in  planes  that  are  parallel 
and  those  that  are  perpendicular  to  the  flow  direction.  In  addition, 
they  use  turbulence  models  which  are  based  on  the  equivalent  sand-grain 
method. 

Flnson  and  his  co-workers  [24,25,26]  have  presented  a  discrete 
element  calculation  method  that  does  not  rely  on  any  sand-grain  concepts. 
The  earlier  work,  Flnson  and  Uu  [24],  does  not  include  the  blockage 
phenomenon.  Their  later  works,  Flnson  and  Clark  [25]  and  Flnson  [26], 
include  both  blockage  factors  and  a  distributed  momentum  sink  term. 
However,  they  do  not  account  for  the  distributed  heat  source/sink  term. 
Their  flow  blockage  factors  are  applied  inconsistently  in  the  equations. 
Also,  their  definition  for  blockage  in  planes  that  are  verpendlcular  to 
the  flow  direction  is  Incorrect.  (This  will  be  discussed  in  Section  III.) 
In  addition,  they  do  not  state  exactly  how  values  for  skin  friction 
coefficient  and  Stanton  number  are  calculated. 

The  work  of  Christoph  and  his  co-workers  [27,28]  follows  the  work 
of  Flnson  very  closely.  They  present  two  methods — an  Integral  method, 
Christoph  [27],  and  a  finite  difference  method,  Christoph  and  Fletcher 
[28].  In  both  papers  they  use  a  mixing-length  model  for  the  turbulent 
shear  stress  that  has  been  modified  to  include  roughness  factors.  In 
the  Integral  method  a  Stanton  numbir  is  calculated  from  the  integrated 
boundary  layer  equation.  Also,  in  the  finite  difference  method  they 
calculate  the  enthalpy  field  over  the  surface.  However,  In  both  papers 
the  Stanton  number  that  is  compared  with  the  experimental  data  Is 
calculated  from  a  correlation  with  the  skin  friction  coefficient.  They 


18 


set  the  heat  transfer  augmentation  due  to  roughness  equal  to  the  square 
root  of  the  skin  friction  augmentation  due  to  roughness. 

Roberson  and  his  co-workers  [29,30,31]  have  worked  continuously  for 
more  chan  twenty  years  with  a  model  for  fully  developed  flow  in 
rough  conduits  that  is  closely  related  to  the  discrete  element  approach. 
Their  model  divides  the  skin  friction  into  viscous  and  element  drag 
components  and  calculates  a  velocity  profile  below  the  crest  of  the 
roughness  elements  based  on  an  assumed  shear  stress  distribution. 

However,  their  model  differs  from  the  discrete  element  model  as  defined 
here  in  that  they  do  not  formulate  a  differential  equation  that  describes 
the  entire  flow  field.  Therefore,  their  insightful  model  is  not  easily 
extended  to  developing  flows  where  boundary  layer  or  Navier-Stokes  based 
calculations  are  riquired. 

The  discrete  element  models  discussed  above,  while  they  are  usually 
insightful,  seem  to  varying  degrees  to  be  patchwork  affairs.  Terms  are 
added  to  the  equations  and  multipliers  supplied  in  an  ad  hoc  manner  to 
satisfy  physical  reasoning.  No  attempt  appears  to  have  been  made  to 
systematically  derive  the  equations  from  first  principles.  In  the 
present  work  a  discrete  element  model  is  derived  for  two-dimensional 
boundary  layer  flow  from  first  principles.  This  model  accounts  for  all 
three  of  the  major  roughness  phenomena  (blockage,  form  drag,  and  distrib¬ 
uted  heat  sink/source) .  The  model  is  calibrated  using  baseline  experi¬ 
mental  data  and  is  compared  to  a  wide  range  of  additional  experimental 
data. 
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SECTION  III 


DERIVATION  OF  THE  DISCRETE  ELEMENT  MODEL 


In  the  discrete  element  approach  the  rough  surface  is  assumed  to  be 
composed  of  discrete  or  distinct  roughness  elements.  Such  a  surface 
could  be  one  composed  of  some  uniform  array  of  cones,  spheres,  hemi¬ 
spheres,  etc.  that  are  glued  or  welded  to  a  smooth  base  surface.  A 
uniform  array  of  elements  with  circular  cross  sections  Is  used  in  this 
section  for  all  derivations  because  It  Is  easily  visualized.  However, 
the  discrete  element  model  that  Is  derived  Is  also  applicable  to  non- 
uniform  arrays  so  long  as  some  adequate  system  exists  to  define  the 
various  parameters  used  to  quantify  the  roughness  geometry.  For  non- 
uniform,  but  descriptive,  arrays  the  geometry  Is  easily  quantified.  For 
example,  consider  an  array  composed  of  rows  of  cones  where  the  distance 
between  cones  In  each  row  and  the  distance  between  rows  is  halved  with 
each  additional  row  until  the  most  dense  packing  Is  obtained.  For 
general  roughness  (l.e.,  where  the  elements  are  random  In  size,  shape, 
and  spacing),  accounting  for  each  Individual  element  Is  no  longer  practi¬ 
cal  and  some  averaging  or  satlstlcal  process  must  be  used.  This  is  not 
a  severe  limitation  on  the  model.  Considerable  spatial  averaging  Is 
already  necessarily  Introduced  by  means  of  the  finite  difference  process. 
In  a  later  section  that  discusses  the  quantification  of  the  roughness 
geometry,  a  system  Is  presented  for  general  roughness. 

In  this  work  only  three-dimensional  or  distributed  (k-type)  rough¬ 
ness  Is  considered.  Graphic  examples  are  the  above  mentioned  arrays  of 
cones,  spheres,  etc.  In  addition,  most  naturally  occurring  rough  sur¬ 
faces  are  of  this  type.  A  rusty  pipe,  a  corroded  turbine  blade,  an 
ablated  heatshleld,  barnacles  on  a  ship's  hull,  a  dull  spray  painted 
surfisce,  and  forests  are  examples  of  distributed  roughness  geometries. 
Two-dimensional  or  strip  (d-type)  roughness  Is  not  considered.  This 
important  type  of  roughness  Is  generally  man-made.  Examples  are  threads 
In  a  pipe,  fiber-resin  composites,  and  rows  of  corn.  The  flow  over  and 
among  strip-type  elements  Is  physically  quite  different  from  that  over 
distributed  roughness  elements  and  requires  a  separate  treatise. 
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In  this  section  the  two-dimensional  time-averaged  boundary  layer 
equations  for  turbulent  flew  ever  a  rough  surface  are  derived.  To  obtain 
the  desired  viquatlons,  a  differential  control  volume  is  imagined  which 
contains  the  pertinent  roughness  effects  determined  from  averaging  over 
several  typical  elements.  Mass,  momentum,  and  energy  balances  are  made 
for  this  control  volume  to  derive  a  system  of  partial  differential 
equations  that  describe  the  flow  over  a  rough  surface. 

As  discussed  in  Section.  1,  empirical  information  must  be  provided 
on  the  nature  of  the  turbulence  and  the  roughness  if  practical  solution 
techniques  are  to  be  used  The  turbulence  model  and  the  roughness  model 
are  described,  and  the  calibration  of  the  roughness  model  is  discussed 
in  detail  in  Section  V. 


1.  DERIVATION  OF  THE  BOUNDARY  LAYER  EQUATIONS 

To  derive  the  boundary  layer  equations  for  thj  conservation  of  mass, 
momentum,  and  energy,  a  differential  control  volume  is  imagined  as  shown 
in  Figure  5.  The  control  volume  is  shown  with  an  exaggerated  length 
in  the  x-dlrectlon  for  easier  visualization  of  the  roughness  effects. 
Using  this  control  volume,  the  appropriate  balances  are  made  to  derive 
the  partial  differential  equations  that  describe  the  boundary-layer  flow 
of  a  fluid  over  a  rough  surface. 

Figure  6  shows  the  control  volume  with  the  rate  of  mass  flow  terms. 
Ferfotnulng  a  steady-state  mass  balance 
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puA  + 
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6x  -  puA 


9puA 
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PvAy  +  6y  -  pvAy 


(6) 


Before  proceeding,  expressions  for  the  area  open  for  fluid  flow  in  the 

x-dlrectlon,  A^,  and  for  the  area  open  for  fluid  flow  in  the  y-dlrectlon, 

A  ,  must  be  developed.  Inspection  of  the  control  volume  reveals  that,  on 

the  average,  not  all  of  either  surface  is  open  for  fluid  flow;  some 

fraction  of  the  surface  is  blocked  by  the  roughness  elements.  This 

phenomenon  is  most  conveniently  expressed  in  terms  of  blockage  factors. 

Let  a  (a  )  be  the  fraction  of  the  surface  that  is  perpendicular  to  the 
X  y 

x-coordinate  (y-coordinate)  that  is  closed  or  blocked  to  flow.  These 

a's  are  hereafter  referred  to  as  the  blockage  factors.  Now  A  and  A 

X  y 
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Figure  6.  Mass 


can  be  expressed  as 


Ajj  -  (1  -  •  1  (7i 

Ay  -  (1  -  ay)6x  •  1  (7, 

where  unit  depth  in  the  z-dlrectlon  Is  assumed. 

Substituting  equations  (7)  into  equation  (6)  and  simplified 
yields  the  desired  continuity  equation 

0  -  [p(l  -  a^)ul  >  Ip-  {p(l  -  ay)v]  (8) 

Figure  7  shows  the  control  volume  with  the  rate  of  momentum  flow 
terms.  Performing  a  steady-state  momentum  balance 
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Again  the  blockage  Is  contained  In  the  terms  A  and  A  as  above.  The  most 

X  y 

important  difference  between  this  balance  and  one  for  flow  over  a  smooth 
surface  is  the  additional  term  F^.  Physically  is  the  force  exerted 
by  the  roughness  elements  on  the  fluid  in  the  control  volume.  This  force 
is  a  result  of  the  shear  and  pressure  distribution  around  the  rou'^hness 
elements  and,  vectorally,  it  is  the  negative  of  the  drag  force  on  the 
elements.  Also,  it  is  convenient  to  chink  of  F^  as  a  distributed  momen¬ 
tum  sink.  That  is,  below  the  crest  of  the  elements  momentum  is  extracted 
volume trlcally  from  the  flow. 

The  term  F^  is  most  conveniently  cast  in  terms  of  a  nondimens lonal 
drag  coefficient,  C^. 

-  i  pC„u^6Ap)  (10) 


where  (SA^)  Is  the  projected  frontal  area  of  the  elements  in  the  control 
volume.  For  elements  of  circular  cross  section 
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Figure  7.  Momentum  Balance  Control  Volume 


6Ap  .  d(y)6y  (  ) 


where  (6x  •  1/ (LZ))  is  the  nuir.ber  of  elements  in  the  control  volume  and 
d(y)  is  the  local  cross  sectional  diameter. 

Substituting  equations  (7),  (10),  and  (11)  into  equation  (9) 
and  simplifying  gives  the  x-momentum  equation 


+  1  pr^2  diil 

2  LZ 


(12a) 


Equation  (12a)  can  be  simplified  by  expanding  the  convective  terms  and 
recalling  equation  (8)  to  give  the  more  usual  form 

p(l  -  a  )u  +  p(l  -  a  )v 
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Figure  8  shows  the  control  volume  with  the  rate  7f  energy  flow 
cenns.  Performing  a  steady-state  energy  balance  with  .‘:he  normal  boundary 
layer  assumptions,  l.e.,  for  shear  work  on  the  x-surface  oW  /dx  ~  0  and 

8  9  X 

for  energy  conducted  through  the  x-surface  3E  ,  /3x  »  0 
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where 
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E  ■  rate  of  energy  convected  across  x-surface  ■ 
cv,x  2 

DuA^(h  -f 

E  ■  rate  of  energy  convected  across  y-surface  * 

y  2 

E  j  ■  rate  of  energy  conducted  across  y-surface  ■  q  A 
cd,y  ^y  y 

W  -  rate  of  shear  work  on  y-surface  ■  TA  u 
8,y  ^  y 

'^drag  "  work  due  to  drag  force  ■  0 

Q  ■  local  heat  transfer  rate  between  element  and  fluid  = 

f  -  '> 


The  two  additional  terms  (over  the  smooth  case)  merit  some  discussion. 
The  work  caused  by  the  drag  (orces  on  the  elements  is  zero  since  the 
velocity  goes  to  zero  at  the  surface  of  an  element.  The  term  Q  is  a  dis¬ 
tributed  heat  sink  (source)  term  .nat  accounts  for  the  volumetric  deple¬ 
tion  (addition)  of  energy  from  (to)  the  flow  because  of  the  presence  of  a 
cold  (hot)  roughness  element.  Here  Q  Is  expressed  in  terms  of  the  local 
Nussclt  nunber,  Nu^,  and  the  roughness  element  temperature,  t^^.  The  term 
Is  the  area  of  the  roughness  elements  in  contact  with  the  fluid  In  the 
control  volume.  For  elements  of  circular  cross  section 

A^  -  Trd(y)6y  (  ^  )  (14) 


Substituting  the  above  along  with  the  expressions  for  A  and  A 

X  y 

(equations  (7a)  and  (7b))  into  the  energy  balance  yields  the  energy 


equation 


•  -  i  1(1  -  0,y)qy)  -  I7  [(1  - 


<'r  -  « 


(15a) 


Equation  (15a)  can  be  expressed  in  the  more  familiar  form  of  static 
enthalpies,  h,  by  combination  with  the  equations  for  continuity  (8) 
momentum  (12) 
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p(l  -  a  )u  +  p(l  -  ot  )v  -1^  ■  -  4“  t(l  -  ot  )q  ] 

X  3x  y  oy  9y  y  y 

-  (1-a  )t|^  +  u-|-  [(1  -  a  )p] 
y  3y  3x  ‘  x 

1  1  j/  \  fTKNu, 

^  2  OV  ^  ^  -U^  ('r  -  '>  (IS'’) 


The  equations  derived  so  far  (continuity,  (8);  x-momentum,  (12); 
and  energy,  (15))  model  steady  laminar  boundary  layer  flow  over  rough 
surfaces.  Before  we  can  practically  apply  these  equations  to  turbulent 
boundary  layers  we  will  have  to  use  the  Reynolds  averaging  process  to 
obtain  the  equations  for  time  mean  turbulent  flow.  First  the  velocities, 
enthalpies,  etc.  are  expressed  as  a  mean  term  plus  a  fluctuating  term 


u  -  U  +  u' 
v  -  V  +  v' 
h  -  H  +  h’ 
p  ■  P  +  p' 
t  -  T  +  t* 


(16a) 

(16b) 

(Ibc) 

(16d) 

(16e) 


Secondly  the  shear  stress,  T,  and  heat  flux,  q  ,  are  expressed  by  the' 
usual  phenomenological  rate  eq  tatlons 


„  3T  K  3h 


(17) 

(18) 


Substituting  equations  (16),  (17),  and  (18)  Into  'he  conservation 
equations  (8),  (12),  and  (15)  and  taking  the  time  average  we 
obtain 


Continuity: 

^  {p(l  -  o.^)Ul  +  Ij  iPd  -  ay)vl  -  0 
x-Momentum : 

(1  -  o.^)pu  If  +  (1  -  a,)PV  If 


(19) 
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(20) 


-  -  [(1  -  a^)Pl  +  I;  [d  -  a^)  (u  f  -  p.'V  )] 


-  7  ft  1  1  77 


,2 


Energy: 


(1  -  a^)pu  f  +  (1  -  a^)pv  |a 


3  r  .  /  K  3H  -TrT\l 

^  L  >' (s^  ). 


+  U  ((1  -  a^)P]  +  (1  -  ay)$ 


,2 


KNu, 


+  TT 


TF  <Tr  - 


(21) 


2  _ 

where  as  usual  Che  terms  3(u'  )/3x  and  3(u'h')/3x  are  neglected.  For  a 
bouxklary  layer  the  mean  dissipation  term,  4>,  Is  approximated  by 


4  a 


3y  '  “  3y 


Inspection  of  equations  (20)  and  (21)  reveals  an  interesting 

2  2 

dependence  on  the  turbulent  Intensity,  u*  /U  ,  in  Che  terms  containing  Che 

drag  coefficient,  C^.  In  general,  Che  turbulent  Intensity  Is  not 

negleglble.  However,  this  dependence  Is  secondary  and  can  be  (to  a  good 

approxloatlon)  Incorporated  Into  the  model  for  the  drag  coefficient,  thus 

— ?  2 

avoiding  a  need  for  another  empirical  closure  for  u'  /U  . 

In  Che  above  we  have  derived  partial  differential  equations  that 
describe  two-dimensional  turbulent  boundary- layer  flow  over  a  rough 
surface.  At  this  point  some  remarks  are  In  order.  The  equations  were 
derived  assuming  distributed  or  three-dimensional  elements.  By  the  very 
nature  of  this  type  of  surface  the  flow  will  be  three-dimensional. 
However,  in  this  model  z-dlrectlon  variations  are  not  considered.  In 
other  words,  all  z-dlrectlon  variations  are  averaged  out.  This  averaging 
out  Is  not  a  severe  limitation  as  long  as  the  control  volume  o.i  which  the 
various  balances  were  made  contains  the  appropriate  amount  of  roughness 
effects  determined  by  averaging  over  several  roughness  elements.  The 
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generally  good  agreement  with  experimental  data  shown  in  later  sections 
reinforces  this  idea. 

The  fact  that  the  control  volume  contains  roughness  effects  deter¬ 
mined  by  averaging  over  several  roughness  elements  is  also  Important  in 
the  flow  direction.  This  averaging  dimension,  in  general,  must  span 
several  elements  if  the  present  model  is  to  make  sense.  This  averaging 
of  the  flow  parameters  over  some  small  but  finite  x-step.  Ax,  is  only  of 
small  consequence,  since  the  most  (and  probably  only)  practical  way  to 
solve  the  partial  differential  equations  is  in  finite  difference  form. 

Before  the  above  equations  can  be  used  to  predict  flow  over  a  rough 
surface  some  system  must  be  developed  to  quantify  the  roughness  geometry. 

2.  QUANTIFICATION  OF  THE  ROUGHNESS  GEOMETRY 

Inspection  of  equations  (19),  (20),  and  (21)  reveals  that  the 

roughness  geometry  is  described  by  the  parameters  a  ,  a  ,  aid  d(y) .  The 

X  y 

parameters  and  Nu^  contain  empirical  information  about  the  drag  and 
heat  transfer  of  the  elements  but  are  considered  here  (as  discussed  in  a 
later  section)  to  be  functions  of  d(y),  or  more  precisely,  the  local 
Reynolds  number,  Re^  ■  U(y)d(y)/v.  First,  we  will  consider  the  case  of  a 
uniform  array  and  then  we  will  consider  a  system  for  general  roughness. 
Here  all  roughness  elements  are  assumed  to  have  cross-sections  that  are 
approximately  circular. 


a.  Uniform  Arrays 

For  uniform  arrays  ';he  cross-sectional  diameter,  d(y) ,  is  the 
same  for  all  of  the  elements  at  a  given  y-location.  Determination  of 
d(y)  is  simply  a  matter  of  considering  the  geometry  of  a  single  element. 
Figure  9  shows  functional  relationships  for  several  common  roughness 
shapes.  Referring  to  Figure  5  we  see  by  inspection  that  is 


a 


y 


(22) 


At  first  glance  i  appears  to  be  a  ■  d/2..  In  fact,  this  is  the 

X  X 

conclusion  of  most  other  workers  [26,28].  However,  this  is  the  maximum 
blockage.  At  most  x-locatlons  the  blockage  is  less,  and  at  some  loca¬ 
tions  the  blockage  may  be  zero.  Here  the  solution  is  to  take  an  average 
x-dlrectlon  blockage.  This  averaging  should  be  natural  in  light  of  the 
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Spherical  Segments 

d(y)  -  2/R^  -  (R  -  k  +  y)^ 


— H 


(d  ny  +  k^  /2k 

o 


Figure  9.  Local  Eleaent  Diameters  for  Cones, 

Spheres  and  Spherical  Segments 

averaging  already  Involved  In  the  finite  difference  nature  of  the  equa¬ 
tions  discussed  above.  For  a  uniform  array  this  averaging  removes  all 
x-dependence  of  a^,  and  the  averaging  length  Is  unimportant  as  long  as 
this  length  spans  several  elements.  For  a  nonunlform  array  some  x- 
dependence  must  be  preserved.  Therefore,  some  averaging  length  must  be 
chosen.  This  length  should  contain  several  elements.  The  natural  choice 
la  a  length  approximately  equal  to  the  x-step  length.  Ax,  used  In  the 
finite  difference  form  of  the  equations. 

The  averaging  process  for  a  uniform  array  of  elements  of  circular 
cross-section  follows.  Since  always  appears  as  (1  -  a the  averaging 
process  will  be  carried  out  on  (1  -  Referring  to  Figure  10 


,  L 

(1  -  a  )  -  T  I  (1  -  a  )dX 

X  ave  L  x 


where 


■  0;  all  X  not  encountering  an  element 

2 

■  — ^^—7 — —  ;  all  X  encountering  an  element 


Substituting  into  equation  (23) 
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Figure  10.  Layout  for  Averaging  a 


(1  -  a  )  ■  T 

X  ave  L 


Completing  the  Integration 


(I  -  d)  ^  jVi  -  j 


dx 


(24) 


(1-a)  •  T  (  I-  '  ^  '  1- 

X  ave  L  \  4*.  / 


Tid^ 

411 


(25) 


It  la  seen  that  (1-a)  -(1-a)  for  any  uniform  array  of  elements 

X  A  V  6  y 

of  Ircular  cross-sections.  This  relationship  also  holds  for  any  array 
of  circular  cross-section  so  long  as  i  and  d(y)  are  not  functions  of  x 
over  the  averaging  distance.  Ax. 

The  reader  may  wonder  why  so  much  time  Is  spent  discussing  the 
simplest  geometry.  There  are  two  strong  reasons  for  considering  uniform 
arrays:  (1)  Insights  can  be  developed  .^y  considering  the  simpler  case 
that  would  be  difficult  to  ^  .certain  from  more  complicated  treatments, 
and  (2)  most  of  the  controlled,  well-defined  experiments  for  flow  over 
rough  surfaces  use  uniform  array  surfaces. 
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b.  General  Roughness 

Before  we  can  quantify  the  geometry  of  a  general  rough  surface, 
we  must  have  some  way  of  collecting  data  on  that  roughness.  The  usual 
way  that  roughness  is  measured  is  by  taking  a  profilometer  trace.  Since 
the  general  roughness  varies  in  both  the  x  and  z  directions,  at  least  two 
traces  are  required  (one  in  the  x  and  one  in  the  z)  .  Figure  11  shows 
typical  traces  and  the  nomenclature  used  below.  Again  we  assume  that  the 
elements  have  circular  cross-sections. 

In  order  to  develop  expressions  for  the  blockage  factors,  we  must 
determine  the  area  blocked  by  the  elements  and  the  associated  plane  areas. 
For  the  xz  plane  the  blockage  factor  is  a.  .  Here  and  in  the  following 
developments  we  will  use  the  z-trace  to  establish  the  element  diameters, 
d(y),  while  the  x-trace  is  used  to  establish  element  density.  Therefore, 


z- Trace 


X- Trace 


Figure  11.  Typical  Profilometer  Traces 


3A 


the  total  xz  area  occupied  by  the  elements  in  the  neighborhood  of  the 
z-trace  is 


A 


tot 


(26) 


where  N^(y)  is  the  number  of  elements  in  the  z-trace  at  a  y-level. 

To  obtain  a  we  need  the  plane  area  associated  with  this  z-trace.  The 
z-dimenslon  of  this  plane  area  is  obviously  L^.  The  x-dimension  for  this 
plane  area  is  the  average  x  spacing  of  Che  elements,  L^/N^(0),  where 
N^(0)  is  the  number  of  elements  at  the  level  y  -  0.  The  xz  blockage  is 
Chen 


a 

y 


TIN  (0) 
X 

4L  L 


X  z 


N^(y) 

I  di^y) 

i«i  ^ 


(27) 


For  the  yz  plane  the  area  occupied  by  an  element  is  where 

d  is  Che  average  diameter  over  some  x  dimension,  I  .  With  a  develop- 
ave  ’  X 

ment  similar  to  chat  which  produced  equation  (25)  is 


d  . 
ave,i 


(28) 


Here  I  - 

X 

blockage 


L^/N^(0)  and  the  associated  plane  area  is  L^6y  giving  the  yz 


a 


X 


TTN  (0) 

X 


4L  L 


X  z 


(29) 


The  form  drag  term  for  general  roughness  is  evaluated  in  a  manner 
similar  to  that  for  uniform  arrays  with  L  replaced  by  L^/N^(0)  and  I 
replaced  by  to  give 


N^(y) 

pu^N  (0)  I  C  d  (y) 

- 2L^ -  ^ 


Likewise,  for  the  local  heat  transfer  between  the  element  and  the  fluid 
is 
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N^(y) 


^KN  (0)  (T_  -  T)  I  Nu 


6x6y  •  1 


For  the  general  roughness  case,  the  boundary  Ir.yer  equations  become 


Continuity: 


[p(l  -  a  )U]  +  -I-  [p(l  -  a  )V]  -  0 


x-Momentum: 


(1  .  |a  *  (1  -  a^)pv  I 


*  If  -  V*  *  ^  [<i  -  v(''  If  ■ 


PU^N  (0) 

"■JiTL 

X  2  1*1 


Energy: 


(1  -  a^)pO  f  +  (1  -  PyjPV  If 

+  U  1^  1(1  -  a^)P]  +  (1  -  tty)? 

pu\  (0) 

"  “2lT 

X  z  !■! 

itra(O)  (I.  -  T)  "z*''* 

*  LL  -  I  "“dl 

X  Z  !■! 

Some  remarks  are  appropriate  here.  The  z  trace  Is  used  to  estab¬ 
lish  roughness  size,  shi|)e,  and  z  spacing.  The  x  trace  Is  used  to  estab- 
the  X  spacing  In  the  neighborhood  of  the  z  trace.  If  significant  x 
variation  of  the  roughness  exists,  several  z  traces  and  associated  x 
traces  would  be  necessary.  Furthermore,  equations  (32)  and  (33) 
reduce  to  equations  (20)  and  (21)  for  uniform  arrays  of  identical 
elements  since  L  ■  L  /N  (0),  i  ■  L  /N  (0)  and  N  (y)  »  N  (0)  for  such  arrays. 

3C  X  Z  Z  Z  Z 
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It  should  be  noted  that  no  effort  has  been  made  to  correct  for  the 
bias  In  profllometer  traces.  This  bias  arises  because  the  trace  has  a 
vanishing  probability  of  passing  through  the  peaks  of  the  roughness 
elements.  Thus,  the  indicated  heights  and  diameters  are  always  on  the 
low  side  of  the  true  values.  No  satisfactory  analysis  of  this  problem 
is  known  to  the  authors. 


3.  EMPIRICAL  CLOSURE 

As  discussed  previously,  empirical  information  must  be  provided  to 
model  the  terms,  u'v'  and  v'h',  that  occur  in  the  Reynolds  averaged 
equations  for  mean  turbulent  flow.  Also,  empirical  information  on  the 
nature  of  the  roughness  effects  must  be  provided  for  a  rough  surface. 


a.  Turbulence  ^k)del 

In  this  work  a  simple  mixing  length  model  with  VanDriest 
damping  is  chosen  to  model  the  Reynolds  stress  term,  u'v'.  That  is. 


-  u'v' 


^  3U 
3y  3y 


y 


iM  1 

T  3y  p 


(34) 


where 

I  »  0.4y(l  -  exp(-y"^/A^)  ] ;  I  <  0.096 
m  m 

I  -  0.096;  I  >  0.096 
m  m 


This  is  exactly  the  lame  model  that  is  so  often  used  for  smooth  surfaces. 
Many  workers  [19,21,32]  attempt  to  model  rough  surfaces  by  adopting  a 
roughness  turbulence  model.  Usually  some  correlation  for  the  damping 
factor.  A***,  as  a  function  of  equivalent  sand-grain  roughness  is  presented. 
Often  some  other  augmentation  of  is  adopted,  for  example,  Christoph 
[27].  The  smooth  wall  turbulence  model  is  adopted  here  for  two  reasons: 
(1)  roughness  effects  are  already  accounted  for  in  the  rough  wall  boundary 
layer  equations  ((19)  -  (21))  and  (2)  as  shown  by  Plmenta  [4]  and 
Coleman  [5]  the  u'v'  and  v't'  terms  are  the  same  nondimenslonally  for 
smooth  and  rough  surfaces. 

To  model  the  term  v'h'  a  constant  turbulent  Prandtl  number,  Pr^  ® 

0.9,  is  assumed.  Then  v'h'  is  given  by 
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(35) 


-  v’h' 


3H  1  -  ^  9H 
Pr^  3y  P  PCp  3y 


b.  Roughness  Model 

The  necessary  empirical  Information  for  rough  surfaces  is 
contained  In  the  models  for  the  drag  coefficient,  and  the  Nusselt 
number,  Nu^.  In  this  work  the  correlations  of  Zukauskas  [33]  for  banks 
of  staggered  cylinders  are  used  as  a  starting  point,  since  at  a  given 
y-locatlon  a  slice  through  an  array  of  elements  with  circular  cross- 
sections  appears  similar  to  a  thin  slice  through  a  bank  of  cylinders. 
Therefore,  the  behavior  of  the  elements  at  a  given  y-locatlon  could  be 
expected  to  exhibit  the  same  general  trends  as  a  bank  of  tubes.  These 
correlations  were  then  modified  by  calibration  with  and  St  data  on 
rough  surfaces.  This  procedure  Is  discussed  In  detail  In 
Section  V  of  this  report. 

The  drag  coefficient  model  developed  is  given  by 


log  Cjj  -  -  0.125  log(Re^)  +  0.375  ,  Re^  <  60,000 


Cp  -  0.6 


,  Re^  >  60,000 


and  the  Nusselt  number  model  by 


«'*d 


2.475  Re. 

a 

1.043  Re^ 
0.963  Re^ 
0.060  Re^ 


0.4p^0.36 

0.5p^0.37 

0.6p^0.36 

0.84p^0.36 


Re^  <  100 

100  <  Re^  <  1000 

1000  <  Re.  <  200,000 
a 

Re^  >  200,000 
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(37) 


It  should  be  noted  that  both  and  Nu^  are  functions  of  y  through 
their  dependence  on  the  local  Reynolds  number,  Re^  ■  Ud(y)/v,  where  both 
U  and  d(y)  are  functions  of  y.  Furthermore,  they  contain  Information  on 
the  element  shape  by  means  of  the  local  diameter,  d(y). 


4.  REMARKS  ON  THE  DISCRETE  ELEMENT  MODEL 

At  this  point  some  remarks  are  In  order.  The  reader  will  recall 
that  in  the  preceding  sections  the  three  regimes  of  flow  over  a  rough 
surface  were  discussed.  No  mention  of  these  regimes  has  been  made  in  the 
development  of  this  model.  All  models  that  use  the  sand-grain  roughness 
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approach  must  take  care  to  distinguish  between  these  regimes  since 
different  models  are  required  for  transitionally  rough  and  fully  rough 
flow.  This  brings  forth  the  added  burden  of  predetermination  of  the 
state  of  the  flow. 

The  present  model  does  not  need  to  make  these  distinctions,  since 
by  its  nature  such  information  is  implicitly  Included.  Therefore,  the 
one  model  applies  to  both  transitionally  rough  and  fully  rough  flows. 

In  Sections  VI  and  VII,  the  present  model  is  compared  with  experimental 
data  for  both  transitionally  rough  and  fully  rough  flows  with  good 
results. 

The  delimitation  of  these  three  regimes  is  usually  made  using  the 

magnitude  of  the  roughness  Reynolds  number,  Re^  .  In  the  present  model 

s 

the  concept  of  sand-grain  roughness,  k^,  has  been  abandoned.  Even 
though  the  concepts  of  fully  rough  and  transitionally  rough  flow  may  be 
artificial  to  some  degree  (witness  Morris'  [2]  discussion  of  an  alternate 
classification),  they  are  a  part  of  our  fluid  mechanics  tradition  that 
(unlike  k  )  we  do  not  yet  wish  to  abandon.  Therefore,  some  new  delimiter 

S 

is  required.  One  candidate  for  this  job  la  the  ratio  of  apnari  shear 
stress  due  to  the  elements  to  total  apparent  shear  stress,  At 

present,  no  limits  are  offered  using  the  new  delimiter  due  to  a  lack  of 
data  in  the  transitionally  rough  regime. 

The  present  discrete  element  model  development  is  for  boundary 
layer  flow.  This  form  was  chosen  for  two  reasons:  (1)  almost  all  of  the 
experimental  data  with  which  the  model  is  compared  are  for  boundary  layer 
flows,  and  C2)  the  relatively  efficient  solution  of  the  boundary  layer 
flow  saved  much  effort  in  the  numerical  experimentation  that  wae  neces¬ 
sary  for  the  callbtac'’x''r  of  the  model.  However,  the  present  discrete 
element  model  should  be  readily  (albeit  tediously)  extendable  to  the 
full  Navler-Stokes  equations. 

5.  DEFINITION  OF  SKIN  FRICTION,  STANTON  NUMBER,  INTEGRAL  LENGTH  SCALES, 

AND  THE  ELEMENT  TEMPERATURE 

At  this  point  formulae  for  the  skin  friction  coefficient  and  Stanton 
number  are  given  in  terms  of  the  discrete  element  model' nomenclature. 

They  are 
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Skin  Friction  Coefficient: 


C 


f 


T 

w  w 


*  i  tl  C 


Stanton  Number: 


^  KiTNu.(IL  -  H) 

-  “y).s,  +  C  — c^rr —  <“>') 


St 


p  U  (H  -  H  ) 
e  e  w  oe 


(38) 


(39) 


The  subscript  w  Indicates  values  evaluated  at  the  wall,  and  e  Indicates 
boundary  layer  edge  conditions.  The  smooth  wall  shear  stress  x  Is  given 

V 

by  -  -  y(3u/ay)|y_Q  and  th^t  heat  flux  is  -  k(3T/8y) 

These  definitions  can  be  formulated  from  physical  reasoning;  howeve*, 
they  also  arise  naturally  In  the  formulation  of  the  integral  boundary 
layer  equations  for  the  discrete  element  model.  These  Integral  equations 
are  der-'ve'’  In  detail  In  Appendix  B  by  Integrating  equations  (19)  - 
(21)  across  ^he  boundary  layer.  During  this  process,  definitions  for 
the  usual  litegral  length  scales  (e.g.,  momentum  thickness,  enthalpy 
thickness,  etc.)  evolve.  The  formulae  for  these  are  given  In 
Appendix  B. 

In  the  development  of  the  energy  equation,  an  element  temperature 
was  Introduced.  This  temperature  In  the  most  general  case  will  have  a 
different  value  than  the  base  wall  temperature,  T^.  The  specification 
of  this  element  temperature  is  rather  Involved  and  is  discussed  In  detail 
In  Appendix  C.  In  the  present  work,  all  of  the  elements  considered  are 
firmly  bonded  and  are  constructed  of  highly  conductive  material.  Therefore, 
here  the  element  temperature  Is  taken  to  be  Identical  to  the  wall 
temperature. 
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SECTION  IV 


METHOD  OF  SOLUTION 


The  discrete  element  equations  derived  In  Section  III  are 
first  transformed  using  a  modified  Illlngsworth  transformation  and  then 
the  transformed  equations  are  solved  using  an  Iterative,  marching. 
Implicit  finite  difference  technique. 


1.  THE  EQUATIONS  AND  BOUNDARY  CONDITIONS 

The  discrete  element  boundary  layer  equations  and  boundary  condi¬ 
tions  as  solved  are: 

Continuity: 


h  (“S/)  + 17  (PS/)  -  0 


Momentum: 


p6u|^+p6v~-p6u  ^ 

X  3x  y  dy  ^e  x  e  dx 


I 


6„(y  +  y») 


3u 

37  I  ^y'^  '  ^T-'  17 
-  T  “S  TI 


Energy : 


+  3y(u  + 


\3y/ 


dU 

p  B  uu  ^ 

ex  e  dx 


1  d  .,3  ,  -  B* 


2  PS  It  "  * 


L£Pr 


The  associated  boundary  conditions  are 
y-0:U-V*0,H-H 

w 

yH-<»;U-^U  ,H-*-H 

•'  e  e 


(40) 


(41) 


(42) 


(43) 

(44) 
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Here  the  nomenclature  has  been  abbreviated  for  compactness.  The 
parametei  3  is  the  fraction  of  the  control  volume  surfaces  open  to  flow, 
that  is,  6^  ■  (1  -  a^)  and  *  (1  -  cXy) .  Though  not  explicitly  indi¬ 
cated,  the  parameters  3,  C^,  Nu^,  and  d  are  all  function  of  y.  In 
addition,  3  is  assumed  not  to  be  a  function  of  x.  For  a  uniform  array  3 
Is  Indeed  only  a  function  of  y,  and  for  most  nonuniform  arrays  3  is  to  a 
good  approximation  only  a  quasi-function  of  x,  i.e.,  d3/dx  =  0.  This 
assumption  is  not  necessary;  however,  it  considerably  simplifies  the 
following  coordinate  transformation.  Also,  all  of  the  data  sets  with 
which  we  compare  in  subsequent  sections  satisfy  this  assumption. 

In  addition  to  the  y-boundary  conditions  stated  above  the  velocity, 
U,  and  the  enthalpy,  H,  profiles  must  be  stated  at  some  initial  x 
position.  This  is  discussed  in  detail  in  Section  IV. 

2.  THE  MODIFIED  ILLINGSWORTH  TRANSFORMATION 

Similar  velocity  and  enthalpy  profiles  exist  for  turbulent  flows 
only  in  Isolated  cases.  However,  considerable  computational  advantages 
occur  if  the  solution  variables  and  coordinate  system  are  chosen  such 
that  similarity  is  approached.  Here  the  Illlngsworth  variables  that 
result  in  similarity  for  the  classical  smooth  wall  laminar  cases  are 
chosen  (for  a  discussion  of  these  cases#  s«e  White  [34]). 

The  following  coordinate  transformation  is  used: 


(45) 

(46) 

(47a) 

(47b) 

(48) 
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Furthermore,  let 


g  -  H/H 


(49) 


From  the  definition  of  the  stream  function,  ip,  continuity,  (40),  is 
satisfied  exactly  and  the  equations  for  conservation  of  momentum  and 
energy  become  (In  the  new  variables)  after  considerable  manipulation 


Momentum: 
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Energy : 


8  +  "'^21*'  ^  ^22®  ^23  ^24  ^  ® 


(51) 
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and 
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The  boundary  conditions  become  In  the  new  variables 
n-0;f'-0,  f-0,g-g^ 

n-^OQ;f'«l,g«l 

The  velocity  components  are  obtained  from  the  expressions 


(52) 

(53) 


(54) 
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(55) 


In  the  above  the  prime  (')  Indicates  the  partial  derivative  with 
respect  to  n,  3(  )/3n. 

At  first  glance  equations  (50),  (51),  and  their  modifiers  appear 
to  be  much  more  intractable  than  their  parents  (equations  (40)  -  (42)). 
However,  we  have  gained  computational  advantages.  We  have  eliminated  one 
variable,  V,  and  we  have  reduced  the  system  of  partial  differential 
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equations  to  a  system  of  quasi-ordinary  differential  equations  with 
quasi-similar  solutions,  f  and  g.  In  other  words,  we  have  sacrificed 
algebraic  compactness  for  computational  considerations. 

3.  THE  IMPLICIT  FINITE  DIFFERENCE  METHOD 

The  development  below  follows  closely  that  for  smooth  surfaces 
referenced  by  Adams  and  Martindale  [35].  Equations  (50)  and  (51) 
are  written  in  the  general  form 

w"  +  Yj^w'  +  Y2W  Y3  +  *  0  (56) 

where  w  *  f  for  the  momentum  equtation  and  w  ■  g  for  the  energy  equation 
with  a  one  to  one  correspondence  for  the  Y's.  That  is,  Y^  “  Y33  for  the 
momentum  equation  and  Y^  “  Y22  t*’®  energy  equation.  Again  the  prime 

indicates  differentiation  with  respect  to  n. 

The  derivatives  in  equation  (56)  are  approximated  by  finite 
difference  formulas  in  the  grid  network  shown  in  Figure  12.  The  grid 

m  m+ 1 


Figure  12.  Layout  of  the  Expanding  Finite 
Difference  Grid 
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spacing  In  Che  ^-direction  Is  a  constant  AC*  A  geometrically  expanding 
grid  Is  used  In  the  n-direction.  This  Is  given  by  the  recursion  formula 

Anijj  -  (57) 

where  R  Is  a  constant  (taken  here  as  R  »  1.063).  The  first  n-step,  Anl^^, 
must  be  specified  and  Is  taken  here  as  An|^  *  0.005.  The  expanding 
n-grld  Is  chosen  because  It  concentrates  the  points  near  the  surface 
where  the  gradients  are  the  steepest.  With  R  ■  1.063,  An|j^  “  0.005  and 
N  ■  120,  43  n  values  are  obtained  with  values  less  than  1  while  the 
largest  n  value  Is  114.  With  this  scheme  on  a  typical  run  about  half  of 
the  n-polnts  are  In  Che  lower  5-lOZ  of  the  boundary  layer. 

The  derivatives  with  respect  to  Che  marching  coordinate,  C,  are 
approximated  by  a  two-point  backward  difference 


rnfl.n 


w 

mfl,n 

AC 


w 

m,n 


(58) 


where  all  values  at  C'Statlon  m  are  known. 

The  derivatives  with  respect  to  Che  normal  coordinate,  n,  are  approxi¬ 
mated  by  using  a  three-point  Lagranglan  interpolation  formula*  The 
dependent  variable  w  Is  approximated  in  the  range  n  i  <  H  <  H  . .  by  a 
Lagranglan  polynomial 


where 
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(59) 


Taking  the  first  and  second  derivatives  of  equation  (59)  and  evaluating 

at  n  *  n  yields  the  finite  difference  formulae 
n 


m+l,n 


2[w  . .  +  Rw  ,  -  (1  +  R)w  ]  , 
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(60) 
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9w 

3n 


m+l,n 


[w  -  R^w  ,  -  (1  +  R)w  ]  . , 
irU _ n-1  _  n  nri-1 

(An)  +  R^(An)  , 

n  n-i 


+  E’ 


(51) 


where 


AnI 


Vi  - 


The  errors  can  be  shown  to  be  [35]  E’  “  -(l/6)w"' (An)  (An)  ,  and  E"  ■ 

n  n— 1 

-(l/12)w""(An)^(An)^^j^.  From  this  we  see  the  usefulness  of  the  expanding 
grid.  Near  the  wall  where  the  derivatives  are  large  we  need  small  values 
of  An  to  give  small  errors.  Far  from  the  wall  the  derivatives  are  small 
and  larger  An's  are  acceptable. 

Substituting  equations  (58),  (60),  and  (61)  into  equation  (56) 
and  rearranging  yields  a  system  of  N-2  algabralc  equations 


where 


A  w  . ,  ,  +  B  w_  . ,  +  C  w  . ,  . , 

n  m+l,n-l  n  nrH,n  n  nrH,n+l 


n-2,  N-1 


The  coefficients  A,  B,  C  and  D  are  given  by 
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where 


(An)^  +  R(An)^  ,  and  S, 
n  n-1  2 


(An)^  +  R"(An) 


n-1* 


The  N-2  equations  are  linear  if  the  terms  A  ,  B  ,  C  ,  D  are  taken 

n  n  n  n 

to  be  constants  at  a  station  n.  Inspection  of  equations  (50)  -  (63) 
reveals  that  these  parameters  are  not  constant.  To  overcome  this 
problem,  the  lollowing  linearization  procedure  is  adopted.  Initial 
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guesses  for  the  velocity  and  enthalpy  profiles  (f  and  g)  are  assumed 
(usually  the  knovm  values  from  the  previous  C~staclon,  m,  are  taken) . 

The  y's  of  equations  (50)  and  (51)  are  evaluated  on  the  basis  of 
these  profiles,  and  the  now  linear  algebraic  equations  (62)  are  solved 
for  the  second  approximation  of  the  profiles.  These  approximations  are 
used  to  obtain  new  approximations  In  an  Iterative  manner  until  the 
difference  between  two  successive  approximations  Is  less  than  a  preset 
percentage. 

This  solution  yields  values  for  f,  not  f.  Inspection  of  equations 

(50)  and  (51)  reveals  that  values  for  f  must  be  obtained  during  the 

iterative  process.  These  values  are  obtained  by  means  of  numerical 

Integration  of  the  f  profile.  In  addition,  values  for  ct  a  I*', 

X  y 

f**',  etc.  are  obtained  by  way  of  the  three-point  Lagrangian  finite 
difference  formulae  discussed  above. 

One  advantage  of  this  formulation  from  a  computational  standpoint 
Is  the  special  tridiagonal  nature  of  the  linearized  equations  (62) . 

These  equations  are  shown  In  matrix  form  in  Figure  13(a).  The  matrix  of 
coefficients  Is  readily  transformed  Into  the  Jordan  canonical  form. 
Figure  13(b),  by  the  relations 
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Figure  13.  Matrix  For*  of  Linearized  Equations:  (a)  Tridiagonal  Form; 
(b)  Canonical  Fom 


where  and  are  known  from  the  boundary  conditions,  equations  (52) 
and  (53).  The  canonical  form  is  easily  solved  for  the  w^'s  by  starting 
at  station  N-1 

"n-I  "  ^-I'^^N-l 


and  back  substituting 


w 

n 


n  -  N-2,  N-3,...2 


(66) 


Here  the  ^-station  subscript,  orfl,  has  been  omitted  for  compactness. 

As  stated  previously,  the  velocity  and  enthalpy  profiles  must  be 
specified  at  some  Initial  ^-station.  One  method  is  to  use  experimentally 
determined  profiles  at  some  location  on  the  particular  body  under  con¬ 
sideration  and  to  then  calculate  the  flow  field  downstream  of  this 
station.  This  s  :heme  has  the  obvious  disadvantage  of  requiring  experi¬ 
mental  data  on  every  surface  that  is  to  be  considered.  Here  the 
following  scheme  is  followed  to  obtain  initial  profiles. 

For  bodies  with  sharp  leading  edges,  such  as  plates,  wedges,  cones, 
etc.,  the  boundary  layer  aquations  are  singular  at  the  leading  edge  and 
meaningful  velocity  and  enthalpy  profiles  cannot  be  obtained.  However, 
starting  profiles  of  the  transformed  variables,  f  and  g,  can  be  obtained. 
The  scheme  is  to  set  the  coefficients,  y'a,  in  equations  (50)  and  (51) 
so  that  the  well-known  Blasius  equations  are  obtained.  These  equations 
are  solved  using  the  implicit  method  discussed  above  with  some  arbitrary 
initial  profile  guess.  For  the  case  of  incompressible,  constant  property, 
laminar  flow  without  pressure  gradients,  this  scheme  yields  the  correct 
profiles  at  every  C-statlon.  For  other  cases  this  scheme  may  introduce 
slight  errors  in  the  solutions  at  the  first  few  C-statlons;  however, 
these  errors  wash  out  after  only  a  few  C-steps. 

For  bodies  with  stagnation  points. 


dU 

U  s  ~ 
e  dx 


x-0 


(67) 


near  x-0.  Equations  (50)  and  (51)  reduce  to  the  classical 
Faulkner-Skan  stagnation  flow  equations  at  x  -  0  (hence  ^  =  0) .  These 
equations  are  then  solved  using  the  implicit  method  with  arbitrary  initial 
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profile  guesses.  The  initial  guesses  used  here  are  f'(r))  ■  1  -  exp(-n) 

and  g(n)  -  g,,  -  (1  -  g  )f'(n). 
w  w 

4.  FULLY  DEVELOPED  FLOWS 


The  above  discussion  Is  concerned  with  Che  more  interesting  case  of 
the  developing  boundary  layer.  However,  irnny  Important  applications  and 
much  of  Che  experimental  data  that  is  used  for  verification  of  Che 
present  model  are  fully  developed  flows.  The  discrete  element  u. 'del  is 
readily  applied  to  such  flows.  For  iully  developed  flow  V  =  0,  and  U 
and  H  are  functions  of  y  only  (i.e.,  J  »  U(y)  and  H  ■  H(y));  therefore, 
the  boundary  layer  equations  (41)  and  (42)  (and,  also,  the  full 
Navier-Stokes  equations)  become  Che  ordinary  differencial  equations 

0  -  -s.  S  *  f  ] 


In  this  form  Che  modified  Illingsworth  transformations  (equations 
(45)-  (49))  are  awkward,  since  there  is  no  previously  known  reference 
velocity  and  enthalpy  (such  as  and  H^) .  However,  no  computational 
advantage  is  gained  by  transforming  the  above  equations,  since  they  are 
already  ordinary  differential  equations.  These  equations  are  readily 
solved  using  the  implicit  finite  difference  method  discussed  above  with 
w  >•  U  and/or  H,  r|~  yi  and  the  appropriate  boundary  conditions. 

5.  COMPUTER  PROGRAM 

The  implicit  finite  difference  method  described  above  has  been 
programed  for  incompressible,  constant  property  boundary  layer  flow  in 
the  Fortran  IV  Plus  language  on  the  Digital  Equipment  Corporation 
PDP-11/34  minicomputer.  Typical  runs  take  5-10  seconds  per  ^-step. 
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The  assumption  of  incompressible,  constant  property  flow  was  used 
here  because  all  of  the  experimental  data  with  which  the  model  is  compared 
in  this  work  meets  this  assumption.  It  should  be  pointed  out  that  the 
model  is  by  no  means  limited  to  incompressible,  constant  property  flows. 

In  fact,  the  discrete  element  model  has  been  coded  and  run  in  the  com¬ 
pressible  variable  property  form  in  a  code  similar  to  the  one  described 
by  Hodge  and  Adams  [21]. 
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SECTION  V 


CALIBRATION  OF  THE  MODEL 


As  with  any  procedure  that  requires  some  empirical  input,  the 
discrete  element  model  must  be  calibrated.  The  calibration  procedure 
here  has  three  parts:  (1)  comparison  of  the  model/computer  program  with 
smooth  wall  results,  (2)  calibration  of  the  momentum  transport  model,  and 
0)  calibration  of  the  energy  transport  model. 

1.  SMOOTH  SURFACE  RUNS 

The  discrete  element  roughness  model  and  computer  program  must 
correctly  predict  the  flow  over  a  smooth  surface  in  the  limit  of  vanishing 
roughness  if  the  model/program  is  to  be  considered  correct.  The  roughness 
program  was  run  for  conditions  corresponding  to  Kearney's  [36]  zero  pres¬ 
sure  gradient  smooth  wall  data  for  U^  >  7  m/s.  To  model  the  smooth 
surface,  the  roughness  model  was  simply  turned  off.i.e.,  the  o's,  d, 
etc.  were  set  to  zero. 

Figure  14  shows  the  comparison  between  the  predictions  and  data  for 
Stanton  number,  St,  and  skin  friction  coefficient,  C^.  The  symbols  are 
the  data  points,  and  the  curves  are  the  predictions.  The  bars  shown  on 
selected  data  points  show  the  reported  experimental  uncertainty  (±10%  for 
C^  and  ±0.0001  for  St).  Inspection  of  the  figure  reveals  that  the 
agreement  is  excellent.  The  predictions  are  everywhere  within  the  data 
uncertainty. 

From  the  above  it  can  then  be  concluded  that  the  model/program 
(exclusive  of  those  parts  dealing  with  the  roughness)  is  correct.  The 
reader  should  note  that  the  present  turbulence  model  has  been  tested  by 
others  for  other  boundary  conditions,  e.g.,  dP/dx  ^  0. 

2.  MOMENTUM  TRANSPORT  MODEL 

The  discrete  element  model  contains  a  term  to  account  for  drag  on 
the  roughness  elements.  This  drag  is  cast  in  terms  of  a  nondlmensional 
drag  coefficient  parameter,  C^.  It  is  this  drag  coefficient  that  con¬ 
tains  the  necessary  empirical  information  on  the  relationship  between 
the  roughness  elements  and  surrounding  flow.  Our  calibration  of  the 
roughness  model  for  momentum  transport  consists  of  the  specification  of 
such  a  relationship. 


53 


0 


0.5 


1.0 


1.5 


2.0 


Figure  1^. 


2.5 


X  [m] 

Comparison  of  Calculations  with  the  Smooth  Wall 
Data  of  Kearney  (36];  *  7  m/s 
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Some  previous  workers,  namely  Flnson  [26]  and  Christoph  and  Fletcher 
[28],  have  used  a  cou^r^nt  value  of  C^.  They  selected  that  value  which 
seemed  to  give  the  best  overall  agreement  with  their  base  data  sets  (i.e., 
Schllchtlng  [3]).  This  scheme  has  a  disadvantage  in  that  must  be 
selected  for  each  roughness  shape.  Finson  and  Clark  [25]  had  good  agree- 
nent  for  surfaces  roughened  with  spheres  and  spherical  segments  using  a 
constant  C^  ■  0.6,  but  they  had  poorer  results  for  the  surface  roughened 
with  cones  when  the  same  value  of  C^^  was  used. 

Lin  and  Bywater  [20]  used  the  correlations  for  banks  of  cylinders 
given  by  Zukauskas  [33].  These  correlations  take  C^  to  be  a  function  of 
the  local  element  Reynolds  number  that  is  calculated  based  on  the  local 
velocity,  U,  and  the  local  element  diameter,  d(y).  However,  they 
evidently  failed  to  realize  that  Zukauskas'  curves  were  expressed  in 
terms  of  a  pressure  drop  parameter  and  not  a  drag  coefficient,  i.e.,  they 
took  the  vertical  coordinate  of  Zukauskas'  Figure  64  on  page  155  of 
reference  [33]  to  be  C^j,  not  the  correct  [(d(y)/f,)Cjj] . 

Here,  also,  the  curves  from  Zukauskaa  [33]  are  used,  but  only  as  a 
starting  point  in  the  specification  of  a  formula  for 

decided  that  to  simplify  the  model,  the  explicit  spacing  (f/d)  dependence 
of  Cq  should  be  dropped.  The  basis  for  this  decision  was  the  fv.ct  that 
a  constant  valus  of  gives  good  agreement  for  a  particular  roughness 
shape,  e.g.,  spheres,  over  a  wide  range  of  roughness  spaclngs.  Further¬ 
more,  spacing  effects  are  already  Included  as  an  integral  part  of  the 
present  discrete  element  model.  Here  and  in  the  section  on  calibration 
of  the  energy  transport  model  that  follows  the  tube  bank  correlations 
are  used  as  a  starting  point  because  in  a  differential  control  volume, 
e.g..  Figure  5,  the  element  slices  appear  to  be  similar  to  a  slice 
through  a  bank  of  cylinders  when  the  elements  have  circular  (or  almost 
circular)  cross-sections. 

To  calibrate  the  model,  various  functions,  C-(Re,),  were  tested  and 

D  a 

the  results  of  the  calculations  were  then  compared  with  a  base  data  set. 
The  function  that  seemed  to  give  the  best  overall  agreement  with  the 
base  data  set  was  then  selected.  Here  the  data  set  of  Schllchtlng  [3] 
was  chosen  as  the  base.  This  data  set  was  chosen  because  it  is  the  most 
comprehensive  data  set  for  flow  over  rough  surfaces  where  the  roughness 
geometry  is  well  defined. 
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As  discussed  previously,  Schlichting  investigated  seven  roughness 

shapes  over  a  wide  range  of  spaclngs  and  over  a  range  of  Reynolds 

numbers  for  a  to*.al  of  124  runs.  In  fact,  this  data  set  has  been  the 

base  for  almost  all  rough  wall  flow  analyses,  including  both  computational 

ones  and  empirical  ones.  Since  in  this  work  we  concern  ourselves  only 

with  three-dimensional  roughness  elements,  only  the  data  for  spheres, 

spherical  segments,  and  cones  are  used  in  the  present  calibration. 

In  the  course  of  this  procedure  a  careful  study  of  Schlichting 's 

report  [3]  revealed  that  his  data  reduction  method  was  flawed.  His  data 

reduction  method  and  the  necessary  corrections  are  discussed  in  detail 

in  Appendix  A.  In  the  present  calibration  of  the  C„(Re.)  model,  the 

D  u 

corrected  skin  friction  coefficient  values  given  in  Appendix  A  are  used 

as  the  data  base. 

Starting  with  functions  (C^  vs.  Re^)  chat  had  the  general  shape  of 
Che  curves  given  by  Zukauskas  [33],  several  different  formulae  were 
tested.  After  some  numerical  experimentation  the  following  function  was 
adopted: 

log  Cjj  -  -0.125  log  (Re^)  +  0.375,  Re^  <  60,000 

Cjj  -  0.6,  Re^  >  60,000  (70) 

Figure  15  shows  a  plot  of  this  function. 

While  equation  (70)  and  Figure  15  indicate  a  relationship  for 
Reynolds  numbers  as  large  as  10^,  the  present  model  has  been  verified 
only  for  roughness  element  Reynolds  numbers  up  to  Re^  ~  25,000.  The 
selection  of  equation  (70)  was  subjective.  That  is,  a  rigorous  optimiza¬ 
tion  procedure  was  not  used.  The  function  chat  seemed  to  give  the  most 
satisfactory  overall  results  for  all  of  Schlichting ' s  cases  was  selected. 

3.  WALL  LOCATION  FOR  MOST  DENSELY  PACKED  SPHERES 

Recalling  the  discussion  in  Section  III,  it  is  seen  Chat  the  discrete 
element  model  divides  the  net  tangential  force  on  the  surface  into  two 
components:  (1)  the  viscous  shear  on  the  smooth  base  wall  and  (2)  the 
drag  on  the  roughness  element.  Most  rough  surfaces  have  fairly  well- 
defined  base  walls.  However,  for  spheres  packed  in  the  most  dense  array 
a  well-defined  base  wall  does  not  exist.  Schlichting  [3]  recognized  this 
problem  in  his  original  article*  arbitrarily  took  the  base  wall  to  be 
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Figure  15.  Roughness  Element  Drag 


at  the  plane  of  the  equators.  The  more  recent  workers  [25,26]  have 
adopted  the  same  approach.  The  general  wisdom  Is  that  the  region  below 
the  equators  offers  such  a  high  flow  resistance  that  the  flow  rate  below 
this  plane  is  negligible.  However,  it  is  more  likely  that  for  densely 
packed  spheres  there  is  a  region  of  recirculating  flow  below  the  plane 
of  the  equators  as  shown  in  Figure  16.  Close  examination  of  the  geometry 
of  the  most  densely  packed  spheres  shown  in  the  figure  rt veals  that  there 
are  most  likely  two  regions  of  recirculating  flow  (one  above  the  equators 
and  one  below) . 


Figure  16.  Schematic  of  Most  Densely  Packed  Spheres 


In  the  absence  of  data  on  this  phenomenon,  the  following  procedure 

was  adopted  to  find  the  correct  (or  at  least  a  good  approximation  thereof) 

effective  wall  location  for  spheres  packed  in  the  most  dense  array. 

Using  the  model  as  calibrated  for  the  preponderance  of  Schllchtlng's  data, 

several  locations  for  the  effective  wall  location  were  tested  and  the  one 

that  gave  the  best  overall  results  for  Schllchtlng's  [3]  skin  frlctldn 

data  for  most  densely  packed  spheres  was  chosen.  The  value  chosen  was 

k  ■  d  /5.  It  is  coincidental  (but  interesting  nonetheless)  that  this 
eii  o 

effective  wall  location  is  almost  identical  to  Schllchtlng's  "melt  down" 
effective  wall  location  discussed  in  Section  II. 
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4.  ENERGY  TRANSPORT  MODEL 

The  discrete  element  energy  transport  model  requires  empirical 

input  in  the  form  of  a  Nusselt  number,  Nu,.  As  with  the  momentum  trans- 

d 

port  model,  the  starting  point  was  the  correlation  reported  by  Zukauskas 
[33]  for  banks  of  tubes.  Lin  and  Bywater  [20]  used  the  same  correlations 
in  tho.i.r  work;  however,  they  used  the  correlations  directly  rather  than 
as  a  starting  point  in  a  calibration  procedure. 

Here  a  model  was  formulated  using  the  correlations  of  Zukauskas  [33] 
as  a  starting  point,  conducting  numerical  experiments  using  modified 
correlations,  and  then  comparing  the  results  of  these  experiments  with 
a  base  data  set.  The  data  set  chosen  as  the  base  was  the  27  m/sec 
experimental  run  by  Pimenta  [4]  at  Stanford  University.  This  data  set 
was  chosen  because  it  is  representative  of  the  very  comprehensive  rough 
surface  heat  transfer  work  carried  out  at  Stanford  [4,5,6,32]. 

The  Stanford  data  were  all  taken  on  a  rough  surface  consisting  of 
spheres  packed  in  the  densest  array.  Therefore,  the  effective  wall 
location  determined  for  Schlichtlng's  most  densely  packed  spheres  was 
used. 

After  the  numerical  experiments,  the  following  Nusselt  number  model 
was  chosen: 


Nu .  - 
d 

2.475 

Re/-Sr0-36  , 

Re^  <  100 

(71a) 

NUd  - 

1.043 

Re,°-V.37  ^ 

100  <  Re^  <  1000 

(71b) 

NUd- 

0.963 

Re  , 

a 

1000  <  Re^  <  200,000 

(71c) 

Nu^  - 

0.060 

Re^0.84pr0.36  ^ 

Re.  >  200,000 
a 

(71d) 

Figure  17  shows  a  plot  of  equations  (71) .  It  should  be  pointed 
out  that  the  model  presented  in  equations  (71)  and  the  figure  have  only 
been  verified  for  Reynolds  numbers  up  to  Re^  ■  1000.  Extrapolations  of 
the  verified  model  using  the  trends  of  Zukauskas'  data  are  presented 
here  for  completeness. 
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SECTION  VI 


FULLY  DEVELOPED  FLOW  RESULTS 


The  discrete  element  model  developed  and  calibrated  in  the  preceding 
sections  has  been  applied  to  two  fully  developed  flow  experiments,  and 
the  results  of  the  calculations  are  compared  below  to  the  results  of  the 
experiments. 

1.  SCHLICHTING'S  EXPERIMENT 

As  discussed  previously,  Schllchtlng  Investigated  several  different 
roughness  shapes  over  a  range  of  spaclngs  and  Reynolds  numbers.  In  this 
section  the  results  of  the  present  model  are  compared  with  the  results  of 
his  experiments.  Since  this  work  Is  concerned  with  three-dimensional 
roughness  elements,  only  the  spheres,  spherical  segments,  and  cones  are 
considered  here. 

As  discussed  In  Section  V  and  In  Appendix  A,  Schllchtlng 's 
[3]  data  reduction  method  for  skin  friction  coefficient  was  flawed,  and 
It  was  necessary  to  correct  this  data.  In  the  following*  all  comparisons 
are  made  with  the  corrected  dt.ta.  It  should  be  noted  that  these  data 
sets  were  used  to  calibrate  the  roughness  model.  The  calculated  results 
cannot,  therefore,  be  considered  as  predictions. 

Table  1  shows  a  comparison  of  the  calculated  skin  friction  coeffi¬ 
cients,  C^,  with  measured  (corrected)  values.  Figures  18  -  20  present 
these  results  In  graphical  form.  The  uncertainty  bands  shown  In  the 
figures  represent  ±10%  about  the  (corrected)  measured  data.  While 
Schllchtlng  did  not  report  data  uncertainties.  It  Is  felt  that  ±10%  Is 
the  smallest  uncertainty  which  can  be  claimed  for  these  data  considering 
the  corrections  which  were  necessary  to  the  originally  reported  values. 

For  the  calculations,  the  discrete  element  model  was  solved  In  the 

proper  channel  coordinates.  The  calculated  values  of  were  normalized 

using  the  measured  value  of  the  maximum  velocity,  U  ,  (l.e.,  - 

2  ypfl  3c  ^ 

2T^/pU^^^  ).  The  maximum  velocity  was  chosen  here  because  Schllchtlng 
did  not  report  average  velocities.  Inspection  of  the  table  (and  figures) 
reveals  good  overall  agreement.  The  calculations  are,  for  the  most  part, 
within  10%  of  the  measured  results  for  11  of  Che  14  plates.  However, 
substantial  disagreement  is  observed  for  three  of  the  plates:  plates  XIII 
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TABLE  1.  COMPARISON  WITH  SCHLICHTING'S  RESULTS 


PLATE 

NO. 

(£/k) 

Re 

X  lO"^ 

U 

max 

(cm/sec) 

MEASURED 

C  X  10^ 

fc 

CALCULATED 

C^  X  10^ 

%  DIFF 

SPHERES 

XII 

100 

321 

6.29 

6.61 

5.1 

(9.75) 

12A 

385 

6.2A 

6.33 

l.A 

162 

A76 

6.A6 

6.09 

-5.7 

190 

5A7 

6.29 

6.00 

-A. 2 

22A 

650 

5.69 

5.90 

3.7 

III 

107 

316 

8.66 

9.61 

10.9 

(A. 88) 

138 

391 

9.18 

9.20 

0.2 

17A 

500 

8.70 

9.22 

6.0 

20A 

568 

9.01 

9.19 

2.0 

251 

70A 

8.95 

9.09 

1.6 

290 

816 

8.99 

8.9A 

-O.A 

I 

107 

310 

18.13 

16.22 

-10.5 

(2.AA) 

132 

38A 

17.87 

15.75 

-11.8 

166 

508 

16. A2 

15. 5A 

-5.5 

195 

566 

16.11 

15.  A3 

-A. 2 

22A 

658 

16.  A2 

15. 3  A 

-6.8 

263 

778 

16.37 

15.32 

-6. A 

II 

lOA 

313 

15.  A6 

15.17 

-1.9 

(1.A6) 

129 

38A 

16.73 

15.60 

-6.7 

166 

500 

17.23 

15. 9A 

-7. A 

186 

586 

17.17 

16.17 

-5.8 

22A 

6A6 

17.57 

16.37 

-6.8 

257 

7A6 

18.00 

16. A9 

-8. A 

V 

98 

311 

8.08 

8.93 

10.5 

(0.87) 

123 

385 

8.AA 

9.02 

6.8 

162 

A98 

8.37 

9.07 

8. A 

190 

585 

8.52 

9.06 

6.3 

21A 

662 

8.56 

9.13 

6.7 

263 

809 

8.86 

9.10 

2.7 

VI 

110 

316 

7.12 

7.85 

10.2 

(A. 86) 

135 

390 

6.82 

7.80 

lA.A 

170 

A91 

7.55 

7.71 

2.1 

200 

566 

7.61 

7.65 

0.5 

23A 

66A 

7.21 

7.62 

5.7 

275 

806 

7.32 

7.63 

A.l 

IV 

107 

325 

13.50 

12.55 

-7.0 

(2. A3) 

129 

391 

12.02 

12.21 

1.6 

186 

572 

12.90 

11.97 

-7.2 

20A 

6A6 

12.59 

11. 8A 

-6.0 

2A5 

751 

12.53 

11.85 

-5. A 
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TABLE  1.  COMPARISON  WITH  SCHLICHTING '  S  RESULTS  (CONCLUDED) 


PLATE  Re  U 

NO-  X  10"^ 

(£/k) _ (cm/sec) 


MEASURED  CALCULATED  %  DIFF 

C-  X  10^  C.  X  10^ 
fc  f 


SPHERICAL  SEGMENTS 


XIII 

115 

314 

4.32 

5.64 

30.5 

(15.4) 

141 

389 

3.72 

5.23 

40.6 

178 

495 

4.37 

5.08 

16.2 

204 

574 

4.40 

4.99 

13.4 

229 

658 

3.70 

4.92 

33.0 

302 

830 

3.33 

4.86 

46.0 

XIV 

135 

380 

4.12 

5.89 

43.0 

(11.5) 

174 

497 

4.57 

5.68 

24.3 

195 

572 

4.87 

5.61 

15.2 

245 

700 

3.68 

5.49 

49.2 

295 

834 

3.50 

5.43 

55.1 

XV 

135 

382 

7.16 

7.70 

7.5 

(7.69) 

170 

502 

7.42 

7.56 

1.9 

190 

564 

7.04 

7.46 

6.0 

234 

687 

6.51 

7.43 

14.1 

288 

817 

6.79 

•  7.40 

9.0 

XIX 

107 

316 

9.93 

10.59 

6.6 

(2.65) 

132 

386  . 

9.75 

10.60 

8.7 

158 

480 

10.27 

10.51 

2.3 

186 

563 

9.08 

10.51 

15.7 

224 

671 

10.21 

10.56 

3.4 

282 

818 

10.25 

10.60 

3.4 

CONES 

XXIII 

117 

321 

4.67 

5.80 

24.0 

(10.7) 

141 

386 

4.89 

5.70 

16.3 

178 

488 

4.69 

5.66 

26.7 

214 

574 

4.84 

5.58 

15.2 

251 

668 

5.12 

5.54 

8.2 

XXIV 

112 

307 

7.81 

7.53 

3.0 

(8.00) 

141 

384 

6.88 

7.18 

4.4 

186 

495 

6.77 

6.87 

1.5 

209 

567 

6.39 

6.77 

5.9 

251 

662 

6.19 

6.65 

7.4 

324 

890 

6.38 

6.54 

2.5 

XXV 

115 

310 

10.32 

9.42 

-8.7 

(5.33) 

145 

388 

8.93 

9.16 

2.6 

178 

476 

8,91 

9.18 

3.0 

209 

564 

9.38 

9.18 

-2.8 

251 

668 

9.34 

9.20 

-1.5 

295 

787 

9.42 

9.18 

-2.5 
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5.0  5.1  5.2  5.3  5.4 

Log  (R«) 

Figure  18.  Comparison  of  Calculations  with  the  Corrected 
Skin  Friction  Data  of  Schlichting  for  Spheres 
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5.0  5.1  5.2  5J  5.4  5.5 

Log(fto) 

Figure  19.  Comparison  of  Calculations  with  the 
Corrected  Skin  Friction  Data  of 
Schllchtlng  for  Spherical  Segments 


^ - 1  f 

5.0  6.1  5.2  5J  5.4  S.5 


Figure  20.  Comparison  of  Calculations  with  the 
Corrected  Skin  Friction  Data  of 
Schllchtlng  for  Cones 
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and  XIV  with  spherical  segments  and  plate  XXIll  with  cones.  For  these 
three  plates  the  calculated  values  of  are  from  8%  to  55%  larger  than 
the  corrected  data.  No  satisfactory  explanation  for  this  disagreement 
has  been  found.  As  discussed  In  the  appendix,  these  three  plates  are 
probably  In  the  transitionally  rough  regime  and  not  the  fully  rough  regime, 
as  originally  reported  by  Schllchtlng.  However,  as  shown  below  In  this 
section  and  In  Section  VII,  the  present  discrete  element  model 
successfully  predicts  other  transitionally  rough  flows.  This  disagree¬ 
ment  is  not  viewed  with  undue  concern,  since  the  model  (as  will  be  shown) 
extrapolr.res  well  to  other  experimental  data  which  were  not  used  In  the 
calibration  of  the  model. 

Figure  21  shows  a  comparison  of  the  calculations  and  data  as  a 
function  of  roughness  density  for  the  spherical  roughness  elements. 
Inspection  of  the  figure  reveals  that  the  discrete  element  model  correctly 
predicts  the  relative  maximum  In  the  — roughness  density  relationship, 

the  maximum  value  of  occurlng  at  a  density  which  Is  somewhat  less  than 
the  maximum  density. 


Figure  21.  The  Effect  of  Roughness  Element  Spacing  on  Skin 
Friction  for  Schllchtlng ' s  Spheres;  Re  =  200,000 
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The  generally  good  agreement  with  Schlichting’s  data  is  not  sur¬ 
prising  since  this  very  same  data  set  was  used  to  calibrate  the  discrete 
element  momentum  transport  model.  The  model  (with  only  changes  that 
account  for  geometric  differences,  e.g.,  from  channel  to  pipe  flow)  has 
been  compared  to  independent  data  sets.  The  pipe  flow  data  of  Chen  [37] 
are  discussed  below,  and  the  extensive  boundary  layer  data  of  the 
Stanford  experiments  are  dl.icussed  in  Section  VII. 

2.  CHEN'S  EXPERIMENT 

Chen  [37]  reported  detailed  turbulence  and  skin  friction  measurements 
for  air  flow  through  a  0.19-meter -diameter  pipe  roughened  with  hemi¬ 
spheres.  He  investigated  three  roughness  densities — l/V.  •  18.5,  10.7,  and 
6.4.  Chen  stated  that  the  first  two  cases  (£./k  ■  18.5  and  10.7)  were  in 
the  transitionally  rough  regime,  and  the  third  (it/k  -  6.4)  was  nearly 

in  the  fully  rough  regime.  The  most  interesting  part  of  Chen's  work 

(from  the  point  of  view  of  the  present  work)  is  the  segregation  of  the 
two  components  of  the  apparent  wall  shear  stress:  (1)  that  due  to  the 
viscous  shear  on  the  smooth  surface  between  the  roughnesi  elements  and 
(2)  that  due  to  the  form  drag  on  the  roughness  elements.  Chen  obtained 

the  form  drag  term  by  measuring  the  force  on  a  single  element  using  a 

force  balance.  The  portion  due  to  the  smooth  surface  was  determined  by 
subtracting  the  roughness  element  drag  component  from  the  total  wall 
shear  stress  which  was  determined  from  pressure  drop  measurements. 

The  discrete  element  model  was  solved  In  the  appropriate  Internal 

circular  coordinates,  and  the  resulting  predictions  were  compared  with 

Chen's  data.  Figures  22  -  24  show  the  comparisons  for  the  skin 

friction  coefficient  and  the  ratio  of  the  smooth  wall  component  to  the 

total  shear  stress.  Here  the  average  velocity  was  used  to  normalize 

2 

the'  skin  friction  coefficient  (C^  *  ^^w^*^^ave  average  velocity 

was  chosen  since  Chen  reported  the  maximum  velocity  only  for  a  limited 
number  of  cases  while  the  average  velocity  was  reported  for  all  cases. 

Inspection  of  the  figures  reveals  very  good  agreement  for  C^  for  the 
first  two  cases  (2./k  ■  18.5  and  10.7).  From  Figures  22  and  23  it  is 
seen  that  the  maximum  disagreement  is  9%.  From  Figure  24,  it  is  seen 
that  fair  agreement  is  obtained  for  the  third  case  (Jl/k  -  6.4);  the 
maximum  disagreement  is  15%.  The  greatest  part  of  this  15%  disagreement 


67 


0.75 


0.50 

0.25 
4.5 

Figure  22. 


5.0 

Log  (Re) 

Comparison  of  Calculations  with  the  Data 
l/V.  «  18.5 


5.5 

Chen; 


68 


4.5  5.0  5.5 


Log  (Re) 

Figure  24.  Comparison  of  Calculations  with  the  Data  of 
Chen;  Jl/k  *  6.4 
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comes  from  a  disagreement  in  average  velocity  and  not  from  the  calculated 
shear  stress,  T^.  In  Figure  24,  it  is  seen  that  the  greatest  disagreement 
occurs  at  the  first  two  points.  Chen  did  report  the  maximum  velocity  for 
the  first  point  (Re  “  60,000).  If  the  maximum  velocity  is  used  to  nor¬ 
malize  both  the  calculated  and  measured  T  ,  the  predicted  value  is  * 

-3  -3  ‘ 

7.31  X  10  and  the  measured  value  is  7.28  x  10 

The  comparisons  of  the  relative  contribution  of  viscous  shear  forces 
between  the  elements  and  the  drag  on  the  elements  are  of  particular 
interest.  One  of  the  major  advantages  of  the  discrete  element  model  is 
that  these  two  forces  and  their  Interactions  are  accounted  for  In  the 
model.  Inspection  of  Figures  22  -  24  reveals  good  agreement  between 
the  predictions  and  data  for  T  /t  .  The  maximum  disagreement  Is  about 
12%  and  the  preponderance  of  the  points  agree  within  5%.  This  agreement 
indicates  that  the  present  discrete  element  model  correctly  Incorporates 
much  of  the  physics  of  the  interaction  between  the  roughness  elements 
and  the  flow. 

Chen  reported  uncertainty  values  of  approximately  ±1%  for  the  data 
presented  in  Figures  22  -  24.  This  small  reported  uncertainty  is 
approximately  represented  by  the  size  of  the  data  point  symbols  in  the 
figures. 
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SECTION  VII 


BOUNDARY-LAYER  FLOW  RESULTS 


The  discrete  element  model  has  been  applied  to  three  of  the  boundary 
layer  experiments  performed  at  Stanford  University.  Two  of  the  experi¬ 
ments  (Healzer  [32]  and  Pimenta  [4])  reported  skin  friction  coefficients 
and  Stanton  numbers  for  both  transitionally  and  fully  rough  zero  pres¬ 
sure  gradient  flow  over  a  constant  temperature  rough  surface.  The  third 
(CoJeman  [5])  included  free  stream  pressure  gradients  and  variable  wall 
lemperature. 

DESCRIPTION  OF  THE  STANFORD  EXPERIMENTS 

The  Stanford  "Roughness  Rig"  consists  of  a  closed  loop  wind  tunnel 
with  air  at  ambient  conditions  as  the  working  fluid.  Free  stream  veloc¬ 
ities  from  9.8  m/sec  to  73.8  m/sec  were  obtained.  This  experimental  set 
up  is  discussed  in  detail  by  Healzer  [32].  The  test  surface  consists  of 
24  0.10-meter  (in  the  flow  direction)  plates  constructed  with  1.27-mm- 
diameter  spheres.  The  surface  was  built  up  of  eleven  tiers  of  spheres 
packed  In  the  most  dense  array.  Each  plate  was  maintained  at  a  constant 
temperature  by  individually  controlled  heating  elements. 

Stanton  numbers  were  obtained  from  an  energy  balance  on  each  of  the 
plates.  The  uncertainty  of  the  Stanton  number  measurements  was  stated 

_3 

to  be  ±0.1  X  10  Stanton  number  units  by  all  three  workers.  One  worker, 
Healzer  [32],  determined  the  skin  friction  coefficient,  C^,  by  differen¬ 
tiating  the  measured  momentum  thickness  distribution;  he  did  not  report 
uncertainty  values.  The  other  workers,  Pimenta  [4]  and  Coleman  [5], 
determined  C^  from  measured  Reynolds  stresses  and  mean  velocity  profiles 
with  a  reported  uncertainty  of  ±10Z. 

One  of  the  concerns  of  the  Stanford  series  of  experiments  was  the 
study  of  the  effects  of  surface  transpiration;  therefore,  both  blown  and 
unblown  runs  were  reported.  Since  transpiration  is  not  considered  in  the 
present  work,  only  the  unblown  cases  are  considered  here. 

2.  APPLICATION  OF  THE  DISCRETE  ELEMENT  MODEL 

The  present  calculations  were  made  using  an  incompressible,  constanf 
property,  discrete  element  boundary  layer  code.  This  more  economical 
code  was  selected  over  the  compressible,  variable  property  code  because 
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the  relatively  low  velocity  and  small  temperature  difference  conditions 
of  the  Stanford  experiments  satisfy  the  Incompressible  and  constant 
property  assumptions.  Furthermore,  the  Stanford  data  reduction  was  made 
based  on  these  assumptions. 

The  roughness  element  drag  coefficient  model  and  Nusselt  number 
model  used  are  those  as  calibrated  In  Section  V.  It  should  be  noted  In 
review  that  the  model  was  calibrated  based  on  Schllchtlng 's  [3] 
channel  flow  data  and  the  Nusselt  number  model  was  calibrated  based  on 
one  of  the  Stanford  runs  (Plmenta  [4],  ■  27  m/sec).  Therefore,  the 

skin  friction  calculations  presented  here  are  true  predictions.  While 
the  Nusselt  number  model  was  calibrated  based  on  one  run  on  the  Stanford 
surface,  its  extrapolation  to  Independent  measurements  under  widely 
varying  conditions  Is  considered  to  be  a  significant  test. 

The  turbulence  model  used  Is  the  mixing  length  and  constant  turbu¬ 
lent  Prandtl  number  model  discussed  In  Section  III,  Some  of  the  Stanford 
runs  were  allowed  to  transition  naturally  from  laminar  to  turbulent  flow. 
For  these  cases  the  transition  point  was  specified  In  the  code,  with 
that  point  selected  based  on  inspection  of  the  Stanton  number  data.  For 
these  cases  a  transition  zone  with  a  length  equal  to  the  laminar  zone 
was  assumed  and  an  Intermlttency  factor  was  employed  In  the  eddy  viscosity 
calculation.  Some  of  the  Stanford  runs  were  tripped  turbulent  by  use  of 
a  boundary  layer  trip.  In  these  cases  turbulent  flow  was  assumed  from 
the  leading  edge.  Virtual  origins  were  not  employed. 

As  stated  previous] y,  the  Stanford  surface  consisted  of  spheres 
arranged  In  the  most  densely  packed  array.  Therefore,  the  effective  base 
wall  location,  as  determined  In  Section  V  using  Schllchtlng 's  data  for 
most  densely  packed  spheres,  was  used. 

In  the  present  calculations,  the  roughness  elements  and  the  surround¬ 
ing  smooth  regions  were  taken  to  be  at  the  same  temperature,  l.e.,  T{^  «  T„. 

3.  CWIPARISONS  WITH  THE  DATA  OF  HEALZER  AND  PIMENTA 

From  the  point  of  view  of  the  present  work,  Healzer's  [32]  and 
Plmenta ’s  [4]  experiments  are  very  similar.  Both  workers  had  a  constant 
free  stream  velocity  and  a  constant  wall  temperature.  Therefore,  they 
are  considered  together.  The  discrete  element  model  has  been  applied  to 
these  experiments,  and  the  skin  friction,  C^,  and  Stanton  number,  St, 
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comparisons  are  discussed  below.  Also,  mean  velocity,  mean  temperature, 
turbulent  shear  stress,- and  turbulent  heat  flux  profile  comparisons  are 
made  for  one  case  (Plmenta's  27  m/sec). 

Figures  25  and  26  show  comparisons  between  the  calculated  and  mea¬ 
sured  values  of  and  St  for  the  two  transitionally  rough  cases 

(Healzer's  U  ■  9.9  m/sec  and  Plmenta's  U  -  15.8  m/sec).  The  bars  on  the 
e  e 

data  points  show  a  ±10Z  span  for  and  the  stated  uncertainty  range  for 
Che  Stantcn  numbers  (±0.001  St  units).  Inspection  of  the  figures  reveals 
that  the  agreement  Is  very  good;  Che  predictions  are  almost  everywhere 
within  the  data  uncertainty. 

Figures  27  -  30  show  comparisons  for  the  fully  rough  case-' 

(Plmenta's  U  ■  27  ra/sec  and  39.6  m/sec  and  Healzer's  U  -57.9  m/sec 
e  e 

and  73.8  m/sec).  Inspection  of  the  figures  reveals  that  the  agreement 
Is  excellent.  Again  the  predictions  are  almost  everywhere  within  the 
experimental  uncertainty. 

Figures  31  and  32  show  comparisons  between  predicted  and  measured 
velocity  and  temperature  profiles.  Inspection  of  the  figures  reveals 
Chat  Che  agreement  Is  very  good.  This  agreement  Is  not  unexpected  since 
such  good  agreement  has  been  demonstrated  for  the  surface  measurables, 
and  Sc.  However,  It  should  be  pointed  out  that  erroneously  calibrated 
discrete  element  models  can  yield  apparently  good  surface  measurable 
agreement  while  at  Che  same  time  giving  poor  boundary  layer  profile 
agreement. 

Figures  33  and  34  show  comparisons  of  measured  and  predicted 
nondlmenslonal  Reynolds  shear  stress  and  turbulent  heat  flux  for  Plmenta's 
27  m/sec  run.  Inspection  of  Che  figures  reveals  that  the  agreement 
Is  excellent.  Such  agreement  bolsters  confidence  In  the  modeling  used 
In  making  the  predictions. 

4.  COMPARISONS  WITH  THE  DATA  OF  COLEMAN 

Coleman  [5]  ran  experiments  or  the  same  apparatus  as  Plmenta  and 
Healzer  with  an  accelerated  free  stream  (favorable  pressure  gradient) 
and/or  variable  wall  temperatures.  His  accelerated  flows  consisted  of 
three  equilibrium  runs  (characterized  by  constant  values  of  Cf  within 
the  equilibrium  zone)  and  one  nonequilibrium  run.  His  variable  wall 
temperature  cases  consisted  of  step  functions,  linear  variations,  end  one 
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Figure  25.  Comparison  of  Calculations  with  the  Data 
of  Healzer;  =  9.75  m/s 
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Figure  26.  Comparison  of  Calculations  with  the  Data  of 

Pimenta;  U  =  15.8  ra/s 
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Figure  27.  Comparison  of  Calculations  with  the  Data  of 
Pimenta;  =  27  m/s 
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Figure  28.  Comparison  of  Calculations  with  the  Data  of 

Pimenta;  U  =  39.6  m/s 
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Figure  29.  Comparison  of  Calculations  with  the  Data  of 
Healzer;  *  58  m/s 
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Figure  30.  Comparison  of  Calculations  with  the  Data  of 
Healzer;  U  =  74  m/s 
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Figure  31.  Typical  Comparison  of  Calculated  and  Measured 
Velocity  Profiles;  Pimenta  »  27  m/s 


Figure  32.  Typical  Comparison  of  Calculated  and  Measured 

Temperature  Profiles;  Pimenta  U  =  27  m/s 
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bi-linear  variation.  Comparisons  with  these  data  are  presented  below. 

Figure  35  shows  the  edge  velocity  variation  for  the  three  equi¬ 
librium  runs.  To  obtain  the  pressure  variation,  dU^/dx,  term  required 
to  solve  the  boundary  layer  equations,  these  data  were  fit  with  piecewise 
linear  functions  using  the  least  squares  method.  Figures  36  -  38 
show  the  comparisons  of  the  skin  friction  and  Stanton  number  predictions 
for  these  flows  with  the  data  for  the  constant  wall  temperature  cases. 
Inspection  of  the  figures  reveals  excellent  agreement.  The  predictions 
are  almost  everywhere  within  the  data  uncertainty.  Figure  39  shows 
the  Stanton  number  comparison  for  the  same  case  as  Figure  35(a)  and  36 
with  step  change  in  wall  temperature  where  the  first  six  plates  (0.61  m) 
were  unheated.  Inspection  of  the  figure  reveals  very  good  agreement. 

Figure  40  shows  the  edge  velocity  variation  for  the  nonequilib¬ 
rium  run.  Again  a  least-square  cur^^e  fit  was  used  to  determine  dlJ^/dx 
in  the  accelerated  region.  Figure  41  shows  the  comparison  of  predic¬ 
tions  and  data  for  this  case  with  a  constant  wall  temperature.  Inspection 
of  the  figure  reveals  that  the  agreement  is  again  excellent.  Figure  42 
shows  Che  Stanton  number  comparison  for  tnls  case  with  Che  first  six 
places  unheaced.  From  this  figure  It  Is  seen  that  Che  agreement  is  good. 

The  last  two  cases  used  for  comparison  are  the  one  with  a  linear 
wall  temperature  variation  and  the  one  with  a  bl-llnear  wall  temperature 
variation.  Figures  43(a)  and  43(b)  show  the  x-variation  of  the  wall 
temperature  for  these  two  cases.  Both  of  these  cases  had  a  constant 
free  stream  velocity  (zero  pressure  gradient).  Figures  44  and  45, 
respectively,  show  the  Stanton  number  comparisons  for  these  two  cases. 
Inspection  of  the  figures  reveals  that  the  agreement  is  excellent. 

Figure  46  shows  a  comparison  of  the  nondimen5ional  shear  stress, 

-  -3 

-u'v',  for  one  of  Coleman's  equilibrium  runs  (K^  =  0.29  x  10  ). 

Inspection  of  the  figuve  reveals  that  the  agreement  Is  reasonable,  with 

the  trend  predicted  co  rectly  but  with  the  predicted  magnitudes  being 

slightly  high. 
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Figure  35.  Edge  Velocity  Distributions  for  Coleman's 
Equilibrium  Pressure  Gradient  Runs 
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Figure  36.  Comparison  of  Calculations  with  the  Data  of 
Coleman;  Kj-  =  0.15  x  10“^;  Constant  Wall 
Temperature 
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Figure  37.  Comparison  of  Calculations  with  the  Data  of 
Coleman;  Kj.  =  0.29  x  10“^;  Constant  Wall 
Temperature 
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Figure  38.  Comparison  of  Calculations  with  the  Data  of 
Coleman;  Kj.  -  0.5  x  IQ-^;  Constant  Wall 
Temperature 
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Figure  39.  Comparison  of  Calculations  with  the  Data  of 
Coleman;  kj.  ■  r.l5  x  10"^;  First  Six  Plates 
Unheated 
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Figure  40.  Edge  Velocity  Distribution  for  Coleman's  Non- 

equilibrium  Pressure  Gradient  Run;  K  -  0.28  x  10"^ 
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Figure  41.  Comparison  of  Calculations  with  the  Data  of 
Coleman;  Nonequilibrium  Run;  Constant  Wall 
Tempeiature 
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Figure  42.  Comparison  of  Calculations  with  the  Data  of 

Coleman;  Nonequilibrium  Run;  First  Six  Plates 
Unheated 
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Figure  43.  Wall  Temperature  Distribution  for  the  Data  of 

Coleman;  Both  Cases  with  Zero  Pressure  Gradient 
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Figure  44.  Comparison  of  Calculations  with  the  Data  of 

Coleman;  Linear  Wall  Temperature  Distribution 
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Figure  ^5.  Comparison  of  Calculations  with  the  Data  of 
Coleman;  Bi-linear  Wall  Temperature 
Distribution 
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Figure  46.  Typical  Comparison  of  Calculated  and  Measured 
Reynolds  Stress  for  Coleman's  Equilibrium 
Pressure  Gradient  Run;  ■  0.29  x  lO”^ 
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SECTION  VIII 


SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

The  present  work  is  sunanarlzed  below.  Also,  particular  conclusions 
and  recommendations  are  enumerated. 

1.  SUMMARY 

The  background  of  predictive  modeling  of  turbulent  flow  over  rough 
surfaces  has  been  discussed.  Of  the  four  broad  predictive  approaches — 
solution  of  the  full,  time-dependent,  Navler-Stokes  equations  in  a 
micro-grid;  solution  of  the  time-averaged  Navler-Stokes  equations  in  a 
micro-grid;  solution  of  the  time-averaged  Navier-Stokes  equations  in  a 
macro-grid;  and  solution  of  the  time-averaged  boundary  layer  equations 
in  a  macro-grid — only  the  last  two  are  practical  with  present  computer 
speed  and  storage  capabilities.  Since  the  macro-grid  is  not  fine  enough 
to  resolve  the  details  of  the  roughness  geometry,  some  form  of  empirical 
input  by  means  of  a  roughness  model  is  required. 

Of  the  two  roughness  sodel  schemes  discussed,  the  equivalent  sand- 
grain  approach  has  several  ^Inadequacies.  It  is  based  totally  on  momentum 
transport  data;  therefore,  there  is  no  physical  reason  for'  the  equiv'ilent 
sand-grain  size  to  correlate  heat  transfer  data.  This  approach  attempts 
to  reduce  a  multi- length  scale  problem  to  a  single  length  scale  problem. 
In  addition,  there  is  no  adequate  means  to  determine  the  sand-grain 
equivalent  for  a  given  rough  surface  without  taking  data  on  that  parti¬ 
cular  surface  or  on  a  similar  surface. 

The  discrete  element  approach  overcomes  many  of  the  inadequacies  of 
the  equivalent  sand-grain  approach.  It  accounts  for  the  momentum  and 
energy  transport  between  the  roughness  elements  and  the  fluid  in  a 
physically  meaningful  manner.  It  accounts  for  the  major  roughness  geome¬ 
try  parameters  (height,  spacing,  and  shape)  and  does  not  rely  on  a  single 
length  scale.  It  Includes  more  basic  physics  of  the  interactions  between 
the  roughness  and  the  flow  and  uses  a  more  general  form  of  empirical 
input. 

In  the  present  work  a  discrete  element  model  for  turbulent  boundary 
layer  flow  over  a  rough  surface  was  rigorously  derived.  The  necessary 
empirical  Information  for  the  roughness  model  was  provided  by  calibrating 
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the  model  using  base  data  sets.  This  calibrated  model  was  then  compared 
with  several  independent  data  sets  covering  both  the  trinsitionally 
rough  and  fully  rough  regimes.  The  results  of  these  comparisons  were,  in 
general,  very  good. 

In  the  course  of  the  present  work  it  was  discovered  that  the  very 
important  data  set  of  Schlichting  [3]  was  flawed  and  that  his  data 
reduction  method  required  correction.  A  discussion  ot  the  reanalysis  of 
this  data  set  and  the  corrected  data  are  presented  in  Appendix  A.  This 
is  a  very  important  consideration  since  Schlichting 's  data  is  the  basis 
for  almost  all  of  the  previous  models  for  flow  over  rough  surfaces. 

2.  CONCLUSIONS 

a.  The  traditional  equivalent  sand-grain  approach  is  restricted 
as  a  predictive  model  for  flow  over  rough  surfaces,  especially  for  heat 
transfer  predictions.  It  may  require  specific  empirical  input  for  a 
given  rough  surface. 

b.  The  discrete  element  model  which  is  derived  in  the  present  work 
for  distributed  three-dimensional  roughness  is  based  on  the  fundamental 
physical  nature  of  the  interactions  between  the  flow  and  a  rough  surface. 
The  resulting  partial  differential  equations  account  for  the  mechanisms 
of  momentum  and  energy  transport  in  a  physically  meaningful  manner. 

c.  In  the  present  model  the  necessary  empirical  information  for 
the  roughness  model  is  of  a  more  general  nature  (than  previous  aiodels) 
and  the  model  can  be  calibrated  for  broad  classes  of  roughnes?;  geometry. 
For  example,  the  calibration  given  in  the  present  work  should  apply  to 
nearly  all  three-dimensional  distributed  roughnesses  whose  cross-sectional 
areas  can  be  reasonably  approximated  by  circles. 

d.  Calculations  using  the  present  discrete  element  model  have  been 
compared  with  118  separate  experimental  runs  on  18  separate  rough  sur¬ 
faces.  The  results  of  these  comparisons  show  that  the  present  discrete 
element  model  is  a  predictive  model  of  acceptable  accuracy  for  both  heat 
transfer  and  skin  friction. 

e.  The  important  data  set  of  Schlichting  [3]  is  flawed.  The 
corrected  values  given  in  Appendix  A  should  be  used  for  any  comparisons 
with  this  data  set.  Any  previous  work  that  is  based  on  the  original 
Schlichting  data  set  should  be  used  cautiously. 
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3.  RECOMMENDED  FUTURE  RESEARCH 

a.  The  experimental  data  base  from  well-controlled  experiments 
with  well-defined  rough  surfaces  is  very  limited.  For  example,  all  of 
the  well-defined  heat  transfer  data  are  based  on  experiments  that  were 
::onducted  on  a  single  rough  surface.  It  is  recommended  that  an  experi- 
rjental  program  be  undertaken  to  establish  a  definitive,  high-quality 
data  base  for  both  uniformly  rough  and  randomly  rough  surfaces. 

b.  The  presence  of  a  rough  surface  affects  not  only  the  value  of 
the  skin  friction  but  also  the  entire  momentum  field  within  the  boundary 
layer.  Therefore,  surface  roughness  could  have  a  strong  effect  on  flow 
separation.  It  is  recommended  that  an  investigation  of  the  effect  of 
surface  roughness  on  flow  separation  be  undertaken. 

c.  The  present  discrete  element  model  is  derived  for  three- 
dlmenslona'  distributed  roughness.  The  important  class  of  two-dimensional 
or  strip-type  roughness  is  characterized  by  different  physical  phenome¬ 
na  and  requires  a  separate  treatise.  It  is  recommended  that  the  discrete 
element  method  be  extended  to  the  two-dimensional  roughness  case. 

d.  Because  of  the  lack  of  well-defined  data  for  higher  speed 
flows,  the  present  model  has  been  verified  only  for  low  speed  flows. 
However,  it  should  be  applicable  to  higher  speed  flows  so  long  as  the 
flow  in  the  immediate  vicinity  of  the  elements  is  subsonic.  It  is 
recomnended  that  calculations  using  the  present  model  be  compared  with 
data  from  experiments  with  high  speed  flows  over  well-defined  rough 
surfaces  as  soon  as  such  data  become  available. 

e.  Usually,  compressible  flow  heat  transfer  data  are  nondlmen- 
slonallzed  using  the  recovery  or  adiabatic  wall  temperature.  The  shear 
rates  in  a  turbulent  boundary  layer  over  a  rough  surface  can  be  signifi¬ 
cantly  different  than  those  for  an  otherwise  similar  flow  over  e  smooth 
surface.  Therefore,  the  adiabatic  wall  temperature  may  be  significantly 
different.  It  is  recommended  that  the  effects  of  surface  roughness  on 
the  adiabatic  wall  temperature  be  investigated. 

f.  Any  profllometer  trace  that  is  used  to  quantify  a  randomly 
rough  surface  will  have  a  built-in  bias,  since  there  is  no  way  to  insure 
that  the  profllometer  stylus  will  encounter  the  crests  of  the  elements. 
Therefore,  it  is  recommended  that  the  statistics  of  random  roughness  be 
investigated. 
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APPENDIX  A 


REEVALUATION  OF  SCHLICHTING’S  SURFACE  ROUGHNESS  EXPERIMENT 

In  Che  course  of  working  with  Che  Schllchclng  [3]  daCa  we  were  led 
CO  quesclon  several  of  Che  assunpclons  used  by  Schllchclng  In  his  daCa 
reducclon.  This  Is  noC  Co  say  ChaC  Schllchclng  made  errors  In  his  analy¬ 
sis  In  1936.  Thac  chls  work  has  been  so  widely  used  by  so  many  researchers 
In  Che  pasC  forCy-plus  years  argues  for  Che  basic  soundness  and  ImporCance 
of  Che  work.  RaCher,  addlC;.onal  daCa  and  analyCical  resulcs  published 
since  1936  now  allow  che  use  of  more  precise  assumpclons  and  enable  a 
reevaluaclon  of  che  original  daca.  The  presenC  avallablllCy  of  che  dlgl- 
cal  compuCer  also  faclllcaces  che  analysis  of  large  daca  secs  and  allows 
more  opclons  Co  be  considered  chan  was  reasonable  wlch  hand  calculaclons . 

In  Chls  appendix,  resulcs  of  che  reevaluaclon  of  Schllchclng 's 
original  daca  for  Che  surfaces  roughened  wlch  spheres,  spherical  segmenCs, 
and  cones  are  presenced.  Ic  Is  shown  ChaC  Che  original  skin  frlcClon 
coefflclenCs  are  higher  Chan  Che  correcced  values  by  amouncs  ranging 
from  0.5%  Co  73%.  Also,  che  equivalent  sand  roughness  values  were 
reevaluated  for  Che  same  seC  of  surfaces.  These  values  are  also  pre¬ 
senced,  and  It  Is  shown  Chac  che  original  values  are  higher  than  the 
corrected  values  by  amounts  ranging  from  26%  co  555%. 

1.  SCHLICHTING'S  EXPERIMENT 

The  experiment  utilized  a  water  tunnel  with  a  test  section  of  40  mm 
height  and  170  mm  width.  A  3.2-m-long  section  wlch  four  smooch  walls 
preceded  che  3.2-m-long  section  in  idtlch  Che  upper  wall  was  rough  and 
Che  side  and  lower  walls  were  smooch.  Fully  developed  channel  flow  thus 
existed  In  the  section  containing  che  rough  wall.  The  14  rough  walls 
wlch  spherical,  spherical  segment,  and  conical  roughness  elements  were 
Chose  of  Interest  In  this  analysis.  The  geometry  and  nomenclature  for 
Che  channel  and  roughness  elements  are  shown  In  Figure  A-1,  and  che 
values  of  che  parameters  for  che  rough  surfaces  are  given  In  Table  A-1. 

The  effective  wall  location  for  a  rough  wall  was  defined  by 
Schllchclng  as  the  location  of  a  "smooth  wall  that  replaces  the  rough 
wall  In  such  a  manner  as  Co  keep  che  fluid  volume  che  same."  Thus,  if  V 
Is  che  volume  of  all  che  roughness  elements  on  a  surface  of  dimensions  L 


103 


J.  J-  AT 


V  'W  V  ^ 

r-f  w  f  f 

1  "  1 

END  VIEW 


Figure  A-1.  Nomenclature  for  Schllchtlng's  Experiment 
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TABLE  A-1.  DESCRIPTION  OF  SCHLICHTING ' S  SURFACE  ROUGHNESS 
AND  CORRECTED  EQUIVALENT  SAND-GRAIN  ROUGHNESS 


Plate 

d 

D 

1 

k 

k’ 

b 

SCHL. 

CORR. 

No. 

(cm) 

(cm) 

(cm) 

(cm) 

(cm) 

(cm) 

k  /k 
s 

k  /k 
s 

SPHERES 

XII 

0.41 

4 

4 

0.41 

a 

3.99 

0.277 

0.120 

III 

0.41 

0 

2 

0.41 

— 

3.99 

0.838 

0.410 

I 

0.41 

1 

0.41 

— 

3.96 

3.07 

2.43 

II 

0.41 

0.6 

0.6 

0.41 

— 

3.88 

3.81 

2.59 

V 

0.41 

0.41 

0.36 

0.41 

— 

3.68 

0.626 

0.378 

VI 

0.21 

1 

1 

0.21 

3.99 

0.819 

0.430 

IV 

0.21 

0.5 

0.5 

0.21 

— 

3.97 

3.61 

2.47 

SPHERICAL  SEGMENTS 

XIII 

0.8 

4 

4 

0.26 

3.99 

0.118 

0.018 

XIV 

0.8 

3 

3 

0.26 

— 

3.99 

0.186 

0.034 

XV 

0.8 

2 

2 

0.26 

— 

3.98 

0.571 

0.278 

XIX 

0.8 

0.8 

0.69 

0.26 

— 

3.85 

1.40 

0.953 

CONES 

XXIII 

0.8 

4 

4 

0.375 

0.425 

3.99 

0.159 

0.046 

XXIV 

0.8 

3 

3 

0.375 

0.425 

3.98 

0.437 

0.122 

XXV 

0.8 

2 

2 

0.375 

0.425 

3.95 

0.996 

0.471 

—  not  applicable 


105 


by  W,  then  Ay  ■  V/LW  is  the  distance  of  the  effective  wall  location  from 
the  smooth  wall  on  which  the  roughness  elements  occur.  This  is  illus¬ 
trated  in  Figure  A-2.  The  mean  channel  height,  b,  was  then  taken  to  be 
the  original  distance  between  smooth  upper  and  lower  walls  (40  mm)  minus 

Ay. 

For  each  of  the  rough  surfaces  used,  either  5  or  6  different 
Reynolds  number  runs  were  made.  Data  taken  were  the  pressure  drop  in  the 
rough  wall  test  section  and  the  velocity  profile  at  the  outlet  plane  of 
the  rough  wall  test  section.  A  pitot  probe  was  used  for  the  velocity 
measurements.  The  flow  rate  was  not  reported,  so  the  reported  maximum 
velocity,  u  ,  is  used  as  the  reference  velocity  in  the  Reynolds  number 
and  skin  friction  coefficient  definitions: 


u  b/v 

(A-1) 

max 

2t/pu^  -  2u*^/u^ 

max  max 

(A-2) 

where  T  is  the  wall  shear  stress  and  u*  is  the  friction  velocity. 

Actually,  for  the  rough  wall  T  is  considered  to  be  the  force  on  the 
wall  In  the  mean  flow  direction  divided  by  the  plan  area  of  the  wall, 
since  the  force  on  the  wall  consists  of  shear  and  form  drag  components. 

In  the  following,  the  words  "wall  shear  stress"  are  used  in  this  sense 
when  referring  to  values  for  the  rough  walls. 

2.  SCHLICHTING'S  EVALUATION  OF 

Schllchtlng  determined  the  ukln  friction  on  the  rough  walls  by  two 
different  methods,  and  then  reported  the  average  of  these  values,  C^.  as 
the  value  of  skin  friction.  The  first  method  he  used  (subscript  1)  was 
based  on  the  pressure  drop  measurements,  while  the  second  method  (sub¬ 
script  2}  was  based  on  the  measured  velocity  profiles  and  a  rough  surface 
law  of  the  wall. 

For  the  first  method,  an  application  of  the  basic  momentum  theorem 
to  the  fully  developed  flow  in  the  channel  yields  the  equation 
Schllchtlng  used 

+  Tg  “  bldP/dx|  (A-3) 

if  one  assumes  the  shear  stresses  on  the  smooth  side  walls  are  negligible. 
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Here,  Che  subscripts  r  and  s  refer  to  the  rough  and  smooth  walls, 
respectively.  This  can  be  recast,  using  the  definition  of  friction 
velocity,  as 


*2  .  *2  b  idPi 

u*  ,  +  u*  ■  —  -T— 
rl  s  p  'dx' 


(A-4) 


The  smooth  wall  friction  velocity  was  determined  by  plotting  the  measured 
velocity  profile  on  the  smooth  wall  in  u  vs.  log  y  coordinates,  deter¬ 
mining  the  slope  n^  graphically,  and  by  comparison  with  the  smooth  wall 
log  lav 

u/u*^  ■  b.5  +  5.75  log  (yu*^/v)  (A-5) 

finding 

u*  -  n  75.75  (A-6) 

s  s 


This  determination,  together  with  the  measured  pressure  drop  data, 
allowed  calculation  of  using  equation  (A-4) .  These  values  are  pre¬ 

sented  in  Table  A-2  in  the  dimensionless  form  for  the  seventy-nine 
runs  made  with  the  rough  surfaces  of  Interest  in  this  study. 

For  the  second  method,  Schllchtlng  used  a  rough  wall  log  law 

u/u*r2  *5.75  log  (y/k)  A  (A-7) 

where  y  is  measured  from  the  effective  wall  location  as  described  pre¬ 
viously  and  k  is  roughness  height.  For  Nikuradse's  fully  rough  sand 

roughness  cases,  A  ■  8.48  and  k  -  k  ,  the  sand  grain  size.  From  the 

s 

velocity  profile  measurements  on  each  rough  wall,  Schllchtlng  knew  u  vs. 
y'  (see  Figure  A-2}.  By  plotting  u  vs.  log  (y*  -  Ay),  determining  the 
slope  n^  graphically,  and  comparing  with  equation  (A-7),  he  was  able  to 
calculate 

u*r2  -  n^/5.75  (A-8) 

These  values  are  presented  in  Table  A-2  in  the  dimensionless  form  C^2* 
Also  shown  in  Table  A-2  are  the  averages  of  and  C^2>  labeled 
Cj  which  were  the  values  reported  by  Schllch  tng. 

3.  REEVALUATION  OF 

As  described  in  the  previous  section,  Schllchtlng  determined  his 
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TABLE  A-2.  COMPARISON  OF  ORIGINAL  AND  CORRECTED 
SKIN  FRICTION  COEFFICIENTS 


SCHLICHTINC  CORRECTED 


Plate 

No. 

(l/k) 

Re 

(xio"^) 

'‘■ax 

(ca/aec) 

'^fl 

(xlO^) 

Cf2 

(xlO^) 

‘^f.av 

(xlO^) 

(xlO^) 

s 

a 

Az 

(cm) 

SPHERES 

XII 

no 

321 

7.60 

7.01 

7.31 

6.29 

75 

0.12 

(9.75) 

124 

385 

7.44 

7.21 

7.33 

6.24 

91 

0.12 

162 

476 

7.52 

7.13 

7.32 

6.46 

114 

0.12 

190 

547 

7.33 

7.31 

7.32 

6.26 

130 

0.12 

224 

650 

6.86 

7.85 

7.35 

5.69 

148 

0.12 

III 

107 

316 

10.22 

11.94 

11.06 

8.66 

299 

0.12 

(4.88) 

138 

391 

10.48 

11.70 

11.08 

9.18 

394 

0.12 

174 

500 

10.06 

12.13 

11.07 

8.70 

486 

0.21 

204 

568 

10.28 

11.83 

11.04 

9.01 

572 

0.21 

251 

704 

10.17 

11.94 

11.04 

8.95 

700 

0.21 

290 

816 

10.19 

11.76 

10.96 

8.99 

806 

0.21 

I 

107 

310 

19.34 

17.16 

18.23 

18.13 

2556 

— 

(2.44) 

132 

384 

18.95 

17.40 

18.17 

17.87 

3118 

0.08 

166 

508 

17.77 

18.70 

18.23 

16.42 

3819 

— 

195 

566 

17.48 

18.94 

18.20 

16.11 

4360 

0.16 

224 

658 

17.69 

18.75 

18.22 

16.42 

5127 

0.16 

263 

778 

17.60 

18.79 

18.19 

16.37 

5944 

0.16 

II 

104 

313 

16.76 

23.80 

20.13 

15.46 

2562 

0.41 

(1.46) 

129 

384 

17.86 

22.92 

20.31 

16.73 

3243 

— 

166 

500 

18.34 

22.07 

20.16 

17.23 

4248 

— 

186 

586 

18.36 

22.01 

20.14 

17.17 

4774 

0.37 

224 

646 

18.70 

21.86 

20.25 

17.57 

5691 

0.41 

257 

746 

19.04 

21  '»3 

20.22 

n.oo 

6708 

— 

V 

98 

311 

9.62 

10.77 

10.19 

8.08 

264 

— 

(0.87) 

123 

385 

9.87 

10.54 

10.20 

8.44 

337 

— 

162 

498 

9.81 

10.64 

10.22 

8.37 

438 

— 

190 

585 

9.95 

10.48 

10.21 

8.52 

525 

— 

214 

662 

9.94 

10.49 

10.21 

8.56 

589 

0.41 

263 

809 

10.17 

10.25 

10.21 

8.86 

739 

0.37 

VI 

no 

316 

8.44 

8.90 

8.67 

7.12 

147 

0.11 

(4.86) 

135 

390 

8.20 

9.18 

8.68 

6.82 

177 

— 

170 

491 

8.70 

8.72 

8.71 

7.55 

234 

— 

200 

566 

3.71 

8.66 

8.69 

7.61 

276 

— 

234 

664 

8.39 

9.02 

8.70 

7.21 

315 

0.11 

275 

806 

8.47 

8.96 

8.71 

7.32 

373 

— 

TV 

107 

325 

14.69 

14.37 

14.53 

13.50 

1154 

— 

(2.43) 

129 

396 

13.44 

15.74 

14.57 

12.02 

1316 

— 

186 

572 

13.96 

15.16 

14.56 

12.90 

1953 

— 

204 

646 

13.68 

15.45 

14.55 

12.59 

2126 

— 

245 

751 

13.62 

15.55 

14.57 

12.53 

2548 

0.13 

—  not  determined. 
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TABLE  A- 2.  COMPARISON  OF  ORIGINAL  AND  CORRECTED 
SKIN  FRICTION  COEFFICIENTS  (CONCLUDED) 


Place 

No. 


Re 


u 

■ax 

(c«/«ec) 


SCHLICHTINC 


CORRECTED 


'“fl 

(xlO^) 


'*£2 

(xlO^) 


f.av 


'•fc 

(xlO^) 


Ax'* 

(c«) 


SPHERICAL  SEGMENTS 


XIII 

115 

314 

5.77 

5.15 

5.46 

4.32 

6 

— 

(15.4) 

141 

389 

5.15 

5.76 

5.43 

3.72 

7 

0.16 

178 

495 

5.56 

5.40 

5.48 

4.37 

10 

— 

204 

574 

5.54 

5.39 

5.46 

4.40 

11 

— 

229 

648 

4.99 

5.91 

5.44 

3.70 

12 

0.16 

302 

830 

4.69 

6.25 

5.44 

3.33 

14 

0.16 

XIV 

135 

380 

5.61 

6.50 

6.05 

4.12 

14 

— 

(11.5) 

174 

497 

5.85 

6.. ’.2 

6.03 

4.57 

18 

0.23 

195 

572 

6.05 

6.08 

6.07 

4.87 

22 

— 

243 

700 

5.10 

7.08 

6.05 

3.68 

23 

— 

295 

834 

4.94 

7.29 

6.06 

3.50 

27 

0.23 

XV 

135 

382 

8.43 

8.70 

8.57 

7.16 

146 

— 

(7.69) 

170 

502 

8.54 

8.53 

8.54 

7.42 

187 

0.18 

190 

564 

8.25 

8.88 

8.56 

7.04 

205 

— 

234 

687 

7.83 

9.37 

8.56 

6.51 

242 

— 

288 

817 

8.02 

9.08 

8.54 

6.79 

202 

0.13 

XIX 

107 

316 

11.29 

11.40 

11.35 

9.93 

484 

— 

(2.65) 

132 

386 

11.13 

11.76 

11.45 

9.75 

590 

— 

158 

480 

11.49 

11.26 

11.38 

10.27 

734 

— 

186 

563 

10.60 

12.23 

11.40 

9.08 

803 

224 

671 

11.44 

11.31 

11.37 

10.21 

1033 

0.26 

282 

818 

11.46 

11.37 

11.42 

10.25 

1296 

— 

CONES 


XXIII 

117 

321 

6.04 

6.83 

6.43 

4.67 

24 

— 

(10.7) 

141 

386 

6.16 

6.72 

6.44 

4.89 

30 

— 

178 

488 

5.98 

6.93 

6.44 

4.69 

37 

0.30 

214 

574 

6.09 

6.83 

6.45 

4.84 

45 

— 

251 

668 

6.28 

6.67 

6.47 

5.1'. 

55 

0.30 

XXIV 

112 

307 

8.60 

8.54 

8.57 

7.31 

79 

-- 

(8.00) 

141 

384 

8.12 

9.12 

R.61 

6.88 

95 

— 

186 

495 

7.94 

9.33 

8.62 

6.77 

123 

— 

209 

567 

7.61 

9.67 

8.61 

6.39 

136 

— 

251 

662 

7.41 

9.87 

8.60 

6.19 

159 

0.38 

324 

890 

7.50 

9.79 

8.61 

6.38 

211 

0.38 

XXV 

115 

310 

11.37 

11.42 

11.40 

10.32 

375 

— 

(5.33) 

145 

388 

10.22 

12.50 

11.33 

8.93 

428 

— 

178 

476 

10.20 

12.62 

11.38 

8.91 

530 

— 

209 

564 

10.54 

12.20 

11.39 

9.38 

638 

— 

251 

668 

10.51 

12.29 

11.39 

9.34 

761 

0.30 

295 

787 

10.55 

12.19 

11.36 

9.42 

909 

0.34 

—  n^t  determined. 
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values  by  neglecting  the  shear  on  the  two  smooth  side  walls.  If  the 
side  wall  shear  Is  Included,  application  of  the  basic  momentum  theorem 
yields 


where  T  Is  the  average  shear  stress  on  the  three  smooth  walls.  The 
s,  av 

value  of  smooth  wall  shear  was  determined  as  described  previously  from  a 
smooth  bottom  wall  velocity  profile  measured  In  a  plane  parallel  to  and 
midway  between  the  side  walls. 

Examination  of  the  data  of  Leutheusser  [38]  for  the  shear  distribu¬ 
tion  around  the  walls  of  a  smooth  channel  of  aspect  ratio  3:1  for  tur¬ 
bulent  fully  developed  flow  shows  that  the  shear  stress  In  this  center 
plane  Is  larger  than  the  average  shear  stress  over  the  channel  perimeter 
by  10  to  20Z,  depending  on  Reynolds  number.  The  present  authors  have 
used  the  following  relationship  between  the  average  shear  stress  and 
T  ,  the  smooth  wall  value  of  the  centerplane  determined  by 

Schllchtlng: 


T 

s,mea8 


1.10 


T 


s,av 


(A-10) 


This  was  chosen  with  the  realization  that,  for  the  larger  aspect  ratio 
used  by  Schllchtlng,  the  Influence  of  the  corner  flows  on  the  average 
shear  stress  should  be  slightly  less  than  In  Leutheusser 's  case.  In 
addition,  of  course,  Leutheusser  had  no  Influence  of  a  rough  upper  wall 
on  his  flows. 

Recasting  (A-9)  and  (A-iO)  In  terms  of  friction  velocities. 


*2  b  dP 

u*  ■  — 
rc  p  dx 


(  W,±.2b  \  ^*2 

I  1.10  W  /  s,meas 


(A-11) 


where  the  rc  subscript  Indicates  a  corrected  rough  wall  value.  All 

values  on  the  right  hand  side  of  (A-11)  were  reported  by  Schllchtlng, 

allowing  calculation  of  for  each  run.  These  values  are  reported  In 

Table  A-2  In  the  dimensionless  form  .  Comparison  of  C,  and  C. 

rc  f,av  rc 

shows  the  originally  reported  values  are  larger  than  the  corrected  ones 
by  0.5Z  to  73Z,  depending  on  Reynolds  number  and  surface  roughness. 

Now  consider  the  velocity  profile  slope  method  of  determining  wall 
shear.  Using  the  z-coordlnate  shown  in  Figure  A-2  to  avoid  confusion 


110 


with  Schlichting's  Ay,  (A-7)  can  be  written  as 

u  ■  5.75  u*p2  ■  Az)  -  5.75  ^  "*■  ^“*^2 

For  a  given  velocity  profile  on  a  given  surface,  the  final  two  terms  in 
(A-12)  are  constants,  so  calculated  from  a  slope  determina¬ 

tion  using  the  measured  data  pairs  (u,z')  if  Az  is  known  a  priori. 

In  using  this  approach,  Schlichting  assumed  that  his  Ay  (as  defined 
previously)  was  equal  to  the  Az  in  (A-12),  which  is  the  wall  shift 
required  to  give  a  velocity  profile  slope  of  5.75  in  u^  vs.  log  z 
coordinates.  Although  the  definition  of  Schlichting's  Ay  is  a  very 
logical  one  when  considered  on  physical  grounds,  it  is  unrelated  to  any 
characteristics  (assumed  or  proven)  of  the  velocity  profile. 

Unfortunately,  for  the  cases  investigated  by  Schlichting,  the  skin 
friction  determined  using  this  technique  is  highly  sensitive  to  the 
value  of  Az.  This  is  illustrated  in  Figure  A-3,  which  shows  the  results 
obtained  for  the  value  of  using  this  technique  and  considering  a 
range  of  Az  values  from  zero  to  the  roughness  element  height.  The 
velocity  profile  data  used  are  for  plate  III,  Re  107,000.  The  slope 
determinations  were  made  using  a  linear  least-squares  regression  in  u 
vs.  log(z'  -  Az)  coordinates.  Plotted  for  comparison  with  the  curve  are 
two  points.  One  is  Schlichting's  (S^^^phical  slope  determination) 
plotted  using  his  value  of  Ay.  The  other  is  the  corrected  value  deter¬ 
mined  from  the  pressure  drop  measurements  (C^^)  plotted  at  the  optimum 
value  of  Az  (discussed  in  the  following  subsection) . 

When  u*  is  not  known  from  an  independent  measurement,  the  question 
of  which  Az  choice  gives  the  best  linear  regression  is  generally  one 
which  cannot  be  answered  unambiguously  when  actual  velocity  profiles  are 
being  evaluated.  This  is  Illustrated  in  Figure  A-4,  which  shows  the 
standard  error  of  estimates  (Schenck,  [39])  vs.  Az/k  for  the  results 
which  are  presented  in  Figure  A-3.  For  this  particular  case,  there  is 
essentially  no  difference  in  the  standard  error  over  the  range  Az/k  from 
0.0  to  0.2,  and  the  difference  for  Az/k  from  0.0  to  0.5  is  only  ±5%  about 
the  average  value  of  standard  error  in  that  range. 

After  considering  the  sensitivity  of  the  velocity  profile  slope 
method  to  the  assumptions  required  to  obtain  a  value  of  u*,  it  was 
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Cue  10* 


Az/k 

Figure  A-3.  Doonscration  of  Influence  of  Effective 

Wall  Location  on  Skin  Friction  Calculated 
From  Velocity  Profile  Slope;  Plate  III, 

Re  -  107,000  ,A-Cf2  .  l-Cf^ 
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Az/k 

Figure  A-4.  Standard  Error  of  Esclmates  for  Velocity 
Profile  Regressions  for  the  Effective 
Wall  Locations  in  Figure  A-3 
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concluded  that  for  the  conditions  of  the  Schlichting  experiment  this 
approach  produces  estimates  of  u*  for  which  the  uncertainties  are  much 
larger  than  for  those  estimates  based  on  pressure  drop  measurements. 
Therefore,  the  corrected  (C^^)  values  based  on  the  pressure  drop  measure¬ 
ments  are  recommended  as  the  best  estimates  of  the  true  values  of  in 
Schlichting ' s  experiment. 

4.  REEVALUATION  OF  EQUIVALENT  SAND  ROUGHNESS 

For  each  of  the  surfaces  he  tested,  Schlichting  determined  an 

equivalent  sand  roughness,  k  .  He  did  this  by  comparing  the  rough  wall 

s 

log  law  in  the  form 

u/u*  -  5.75  log(y/k)  +  A  (A-13) 

to  the  form  Nllcuradse  reported  for  his  velocity  profile  data  in  the 
fully  rough  regime 

u/u*  ■  5.75  log(y/k  )  +  8.48  (A-14) 

s 

tt  should  be  noted  that  in  (A-13), 

y  -  y’  -  Ay  (A-15) 

while  Nikuradse  [1]  did  not  explicitly  define  the  origin  of  his  y  coordi¬ 
nate  used  in  (A-14).  Schlichting  set  (A-13)  and  (A-14)  equal,  assumed 
all  of  his  data  were  in  the  fully  rough  regime,  and  obtained 

5.75  log(k  /k)  -  8.48  -  A  (A-16) 

8 

By  computing  an  average  value  of  A  for  each  velocity  profile  using 

(A-13)  and  then  calculating  a  mean  value  of  A  from  all  the  profiles  on 

each  plate,  he  was  then  able  to  use  (A-16)  to  solve  for  a  k  /k  value  for 

s 

each  of  the  rough  surfaces  he  tested.  These  values  are  listed  in  Table 
A-1  for  the  surfaces  of  Interest  in  the  present  study. 

In  view  of  the  discussion  in  the  preceding  sections,  the  use  of 
u*  and  Ay  in  (A-13)  are  open  to  serious  question.  It  is  more  appro- 

IT  p  8V 

prlate  to  use  the  corrected  friction  velocity  and  the  wall  shift  Az  such 
that 

u/u*^^  ■  5.75  log[(z'  -  Az)/k]  +  A  (A-17) 
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If  (A-17)  Is  used  in  determining  an  equivalent  sand  roughness,  then  Az 
values  must  be  determined  before  the  computation  can  proceed. 

Monln  and  Yaglom  [9]  show  chat,  within  the  logarithmic  layer,  Che 
quantity 

z  ■  (z'  -  Az)  exp(-Ku/u*)  (A-18) 

0 

is  Independent  of  distance  from  the  wall.  Thus,  if  the  friction  velocity 
and  Che  data  pairs  (u,z')  in  a  velocity  profile  are  known,  z^  can  be 
determined  as  a  function  of  (z*  -  Az)  for  various  Az  values.  The  opti¬ 
mum  Az  value  is  Chen  the  one  which  gives  values  of  z^  which  are  the 
closest  to  being  constant  with  (z*  -  Az) .  The  optimum  Az  values  were 
determined  in  this  study  using  Che  above  approach  and  a  criterion  of 
minimum  standard  error  of  estimates  from  a  lines'*  least  squares  regres¬ 
sion  of  z^  on  (z*  -  Az)  with  zero  slope.  The  Karman  constant,  <,  was 
taken  as  O.AO  for  consistency  with  the  slope  in  (A-17). 

In  using  this  procedure,  the  velocity  profile  points  must  be  read 
from  Schllchtlng's  [3]  figures.  The  symbols  for  the  different  profiles 
are  not  distinguishable  in  most  cases  in  the  NACA  English  translation, 
so  the  original  German  version  was  used.  Shown  in  Table  A-2  are  the 
calculated  "optimum"  Az  values  for  each  of  the  velocity  profiles  which 
could  be  read  from  the  figures. 

For  each  of  the  profiles  for  which  Az  was  calculated,  a  linear 

least  squares  regression  of  Che  form  of  (A-17)  was  used  Co  determine  a 

value  of  A.  Then  (A-16)  was  used  to  calculate  the  corresponding  k  /k 

s 

value.  The  mean  values  of  k  /k  (shown  in  Table  A-2)  were  then  deter- 

s 

mined  for  each  surface.  Comparisons  show  that  Schllchtlng's  original 
k^/k  values  are  higher  than  the  correct  ed  values  by  26Z  to  5S5Z. 

In  using  (A-16)  to  determine  It  is  assumed  that  (z*  -  Az)  for 

Schllchtlng's  data  and  Nikuradse's  y  are  equivalent.  In  other  words,  the 
origins  for  the  coordinate  systems  must  be  those  which  give  a  slope  of 
3.75  in  Che  log  regions  for  both  sets  of  data.  It  has  been  shown  pre- 
vlously  Chat  Az  meets  this  requirement  by  its  definition.  However,  since 
Nlkuradse  did  not  explicitly  define  his  origin  for  y,  Che  sensitivity  of 
his  results  to  a  wall  shift  must  be  considered. 

Fortunately,  Nikuradse's  results  are  much  less  sensitive  to  assumed 
wall  position  than  are  those  of  Schllchtlng.  Nlkuradse  tested 
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configurations  over  a  range  of  r/k  from  15  to  507,  where  r  is  pipe  radius 
and  k  is  the  average  sand  size.  The  analogous  parameter  In  Schlichting ' s 
configurations,  b2/k,  ranged  from  6.1  to  13.8  for  the  surfaces  considered 
in  this  study.  Schlichting  defined  62  as  the  distance  from  the  rough 
wall  to  the  maximum  in  the  velocity  profile  across  the  channel.  If  it 
is  assumed  that  the  virtual  or  effective  wall  location  lies  between  the 
bottom  and  the  top  of  the  roughness  elements,  then  the  points  in 
Nlkuradse's  velocity  profiles  are  much  less  effected  than  are  Schlichting 's, 
since  as  z*  becomes  much  greater  than  Az,  log(z'  -  Az)  approaches  log(z'). 

Using  the  corrected  k^  values  for  each  surface  and  the  corrected 
friction  velocities  a  roughness  Reynolds  number 

Rej^  -  (A-19) 

was  calculated  for  each  of  Schlichting 's  runs  which  were  considered. 

These  values  are  recorded  in  Table  A-2.  It  should  be  noted  that  for  all 
runs  using  plates  XIII,  XIV,  and  XXIII,  the  calculated  roughness  Reynolds 
numbers  are  in  the  transitionally  rough  regime,  not  in  the  fully  rough 
regime  as  assumed  by  Schlichting.  Of  course,  the  method  by  which  the  k 

s 

values  were  determined  assumes  fully  rough  flow,  so  the  meaning  of  the 
k^  values  for  these  three  plates  is  not  the  same  as  for  those  plates 
on  which  fully  rough  flow  existed. 
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APPENDIX  B 


INTEGRAL  BOUNDARY  LAYER  EQUATIONS 
FOR  FLOW  OVER  A  ROUGH  WALL 


The  Integral  forms  of  the  boundary  layer  equations  are  often  very 
useful  when  approximate  boundary  layer  analyses  are  desired.  Also,  their 
development  provides  natural  definitions  for  skin  friction  coefficient, 
Stanton  number,  and  the  usual  Integral  length  scales  such  as  momentum 
thickness,  displacement  thickness,  enthalpy  thickness,  etc. 

For  boundary  layer  flow  over  a  rough  surfsce,  the  following  equations 
have  been  derived  In  the  text.  Furthermore,  we  recall  that  for  most 


surfaces  of  Interest  (1  -  a  ) 

X  ave 

Continuity: 


(1  -  Uy)  -  (1  -  a). 


1^  [p(l  -  a)u]  +  tp(l  -  a)v] 


(B-1) 


Momentum: 


p(l  -  “)u  1^  +  P(1  -  a)v  1^  -  -  ly  [(1  -  a)T] 
+  (1  -  a)p^U^  ^  i  pc^u^ 


(B-2) 


Energy: 


p(l  -  a)u  1^  +  p(l  -  a)v  1^  -  -  1^  [(1  -  a)q] 
a  dU 

-  (1  -  f  -  “Pe<"  -  “^^e  ^ 

1  ^  3d  . 

2  PS'*  n  — mi — 


(B-3) 


First,  we  will  consider  the  integral  momentum  equation.  Multiplying 
continuity  by  (u  -  U^) ,  adding  it  to  the  momentum  equation  and  rearranging 
yields 

a  dU 

-  [p(l  -  a)u(U^  -  u)l  -  (1  -  a)(p^U^  -  pu) 


+  [p(l  -  a)uv)  -  Ug 


9p(l  -  a)v 
3y 
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(B-A) 


((l-o.)Tl  -ipc/fj 


Integrating  equation  (B-A)  across  the  boundary  layer  with  respect  to  y 
and  recalling  that  u^  -  ■  0,  ■  -x^  _  ^  and  *  0  yields 


-  I-  [u 

dx  c 


r  p(i  -  a) 


dU 

djT  ^e 


“e  (  '  ■  U 


+  U,  (1  -  a)^„v„ 


-  -  (1  -  i  C 


(B-5) 


If  we  make  the  following  definitions 


(B-6) 


momentum  thickness  -  9  «  /*  (1  -  a)  (  1  -  ■jj-  ]  dy 

o  ^e^e  \  ^e  / 

displacement  thickness  •  6*  ■  /*  (1  -  a)  |  1  -  \  dy  (B-7) 

o  V  ^e  e  / 


skin  friction  coefficient  “  y  " 


equation  (B-5)  takes  the  familiar  form 


P  U 
e  e 


^  +  (20  +  6*)  ^  ^ 
dx  U  dx  p  dx 

e  e 

p  (1  -  a  )v 
w _ w  w 

P  U 
e  e 


(B-9) 


In  the  last  term  of  equation  (B-9),  the  factor  (1  -  Is  the 

average  transpiration  velocity  based  on  the  total  plan  area. 

Next  we  will  consider  the  energy  equation.  We  multiply  momentum, 

2 

(B-2),  by  u;  define  stagnation  enthalpy,  h^  ■  h  +  u  /2;  multiply  continuity 
by  h^;  and  combine  these  two  equations  with  energy  to  obtain 
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(B-10) 


1^  [p(l  -  a)uh^]  +  (p(l  -  a)vh^] 

-  -  1^  [(1  -  a)ql  -  1^  t(l  -  a)uT] 
TryNu^(hj^  -  h) 


Integrating  equation  (B-10)  across  the  boundary  layer  with  respect  to  y 
and  recalling  u  ■  ■  q^  ■  0  we  get 


(1  -  + r 


TryNu^(h^  -  h) 


-  r  T"  -  h  ) 

j  33j  '  o  oe' 


-  p(l  -  a)  V  (h  -  h  ) 
w  w  o  oe  w 


(B-11) 


where  we  have,  without  any  loss  In  generality,  made  the  stagnation 
enthalpy  at  the  boundary  layer  edge,  h^^,  the  datum  for  stagnation 
enthalpy.  Normalizing  equation  (B-11)  with  expanding 

the  Integral  we  get 


<1  -  *  r 


TryNu^(hj^  -  h) 


e  e  w  e 


p  (1-a)  V  (h  -h  ) 
^w  w  w  0  oe  w 

p  U  (h  -h  ) 
e  e  w  e 


1  d  ,»  . P(l-a)u(h-h^) 

'*^pU(h  -  h)dx^  ^e^e^^w  ^e^  p  U  (h  -h  ) 

~  A  '  ».«  A*  /%  '  A  A  '  r.t  A 


e  e  w  e 


e  e  w  e 


1  1  d  3  p(l-o)u(U  ^-u^) 

-  T  p  U  (i  -h  )  S  /  »e".  - —3^ - -*>' 

e  e  w  e  o  p  U 

e  e 

If  we  make  the  following  definitions 

^  imNu^(h^  -  M  ^ 

Stanton  number  ■  St  ■  ”  ^  -  (B-13) 

"e  e'  w 


h  -  h 


enthalpy  thickness  “  \  (l~ot)  ^  dy  (B-14) 

^e  e  w  "e 
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dissipation  thickness  ■  *5^  ■  /”  (1-a)  [  1  ~  |  ‘^y  (B-15) 


e  e 


Eckert  number  -  E 


c  2(h  -  h  ) 
w  e 


(B-16) 


equation  (B-12)  takes  the  familiar  form 

[1  dp  ,  dU  , 

e  e  w 


d(h  -  h  ) 
w  e 

dx 


] 


d6.  d6 .  r  E  dp  3E  dU  ■ 

_ h  _  „  _ d  _  p  _c  _ e  _ c  _ e 

dx  c  dx  d  p  dx  U  dx 

L  e  e  J 


P  V  p 

-  (1  -  a)  +  (1  -  a)  V  ^  E 

w  p  U  w  w  p  c 

e  e  e 


(B-17) 


As  an  alternate  formulation,  we  can  normalize  equation  (B-11)  with 

reapect  to  the  stagnation  enthalpy  defect,  p  U  (h  -  h  ),  and  define  a 

e  e  ow  oe 

stagnation  Stanton  number  and  enthalpy  thickness 


S^o- 


Truiiu  (IL  -  h) 

(1  -  +  r - pto —  dy 


p  U  (h  -  h  ) 
'^e  e  ow  oe 


(B-18) 


V  ■  r  (1  ^ 


h  -  h 
o  oe 

p  U  h  -  h 
e  e  ow  oe 


dy 


(B-19) 


to  get  the  other  familiar  form 
d6^ 


St  -  — ; — 
0  dx 


+  6. 


,  dp  ,  dU 
1  e  1 _ e 

p  dx  U  dx 

e  e 


d(h  -  h  ) 
ow  oe 


h  -  h 
ow  oe 


dx 


-  (1  -  a) 


P  V 

w  w 

w  p  U 
e  e 


(B-20) 
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APPENDIX  C 


ON  THE  PROBLEM  OF  ROUGHNESS  ELEMENT  TEMPERATURE 

In  the  derivation  of  the  energy  equation  for  the  rough  surface  bound¬ 
ary  layer  It  was  assumed  that  the  temperature  of  the  roughness  element  was 
known.  The  problem  of  defining  this  temperature  Is  addressed  below. 

The  temperature  distribution  within  the  element  will.  In  general,  be 
three-dimensional.  This  temperature  distribution  could  be  found  by 
solving  the  three-dimensional  conduction  equation  with  appropriate  boundary 
conditions.  This  would  be  a  significant  undertaking  on  the  order  of 
solving  the  two-dimensional  boundary  layer  equations.  Also,  since  the 
temperature  profile  must  be  determined  at  each  Iteration  of  the  boundary 
layer  equations,  the  required  computations  would  be  Increased  by  an  order 
of  magnitude. 

Fortunately,  the  exact  temperature  distribution  In  the  element  Is  of 
little  concern,  since  the  major  concern  la  the  total  heat  transfer  rate 
between  the  fluid  and  the  substrate.  This  heat  transfer  is  affected  by 
the  roughness  elements  In  two  ways: 

(1)  The  elements  Influence  the  nature  of  the  turbulence  and  the 
momentum  field,  thereby  influencing  the  energy  field  In  the 
boundary  layer. 

(2)  The  elements  Influence  the  path  for  conduction  for  the  fluid 
to  the  substrate  In  such  a  way  as  to: 

(a)  Increase  the  effective  surface  area  (fin 
effect);  or 

(b)  Insulate  the  surface. 

Effec."  (1)  must  be  accounted  for  In  the  models  of  the  turbulence  and  the 
momentum  field.  Effect  (2)  Is  Included  In  the  model  of  the  energy  field 
by  way  of  the  blockage  terms  and  the  distributed  heat  source  (sink)  terms. 
It  Is  the  heat  source  (sink)  term  that  requires  a  definition  of  roughness 
element  temperatures. 

The  lumped  mass  approach  described  below  gives  a  simple  method  to 
define  an  element  temperature  that  accounts  for  the  Important  effects  of 
conduction  through  the  element  and  contact  resistance  at  the  boundary 
between  the  element  and  the  surface. 
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Figure  C-1  shows  ?.  schenatlc  diagram  of  the  element  and  the  nomencla 
ture  used  In  this  appendix.  It  Is  assumed  that  the  mass  of  the  element 
Is  lumped  at  the  centroid.  For  steady  state 


q 


I] 


(C-1) 


Solving  for  T_  yields 

K 

Ty  +  +  R2)  /q  (<i)HT(y)dy 

^  1  +  7r(Rj^  +  R2)  (d)Hdy 

o 

where 


and 


0 


/q  (d^)ydy 

'.i^dy 


(C-2) 


(C-3) 


(C-4) 


The  reader  should  note  that  In  the  limit,  as  R.  +  R.  0,  T.  T  ; 
80,  for  highly  conductive,  firmly  bound  elements  the  element  temperature 
approaches  the  wall  temperature.  On  the  other  hand,  as  R^  +  R^  l.e. 

weakly  conductive  and/or  loosely  bound  elements,  the  element  temperature 
approaches  a  weighted  average  of  the  surrounding  fluid  temperature 


as  Rj^  +  R^  -*■  ®, 


/J  (d)HT(y)dy 
(d)Hdy 


(C-5) 
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q  ■  rate  of  heat  transfer  between  element  and  fluid 

■  conductive  resistance  between  centroid  and  surface 
R2  *■  contact  resistance  at  boundary 

k  >  element  height 

■  "temperature"  of  the  element 
T  *  temperature  of  the  substrate 

V 

*  T(y)  >  temperature  of  the  fluid 

y  "  location  of  the  centroid 
c  « 

H  convective  heat  transfer  coefficient 
d  -  local  element  diameter 

Figure  C-1.  Roughness  Element  Heat  Transfer  Schematic  and  Nomenclature 
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